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ABSTRACT 


The  use  of  a  self  adaptive  technique  to  compensate  for  the  varying 
parameters  of  a  flight  control  system  is  desirable  for  high  performance 
aircraft  and  space  vehicles.  The  technique  discussed  in  this  study  con¬ 
sists  of  using  a  sinusoidal  signal  to  measure  the  change  in  the  parameters 
and  of  compensating  for  the  change  by  a  gain  adjustment.  Analog  and  digital 
computer  methods  are  used  to  investigate  the  performance  of  the  system  and 
the  accuracy  of  the  Amplitude  Closed  Loop  Response  Criterion  for  measuring 
damping  ratio.  The  effects  of  varying  airframe  flexibility,  center  of 
gravity,  frequency  of  the  sinusoidal  signal*  and  the  Reference  signal  value 
are  included.  The  applicability  of  the  self  adaptive  technique  for  a 
particular  vehicle  depends  on  t he  dynamic  characteristics  of  the  airframe 
and  on  the  desired  rra  ximum  variations  of  the  performance  characteristics. 
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CHAPTER  I 


INTRODUCTION 

The  dynamic  characteristics  of  an  airframe  vary  when  changes  in 
speed,  altitude,  position  of  center  or  gravity,  and  dynamic  pressure 
occur.  The  variation  in  the  dynamic  characteristics  cause  variations 
in  the  response  of  the  airframe  to  commands  from  a  pilot  or  program  un¬ 
less  some  method  of  compensation  is  used.  Many  autopilots  employ  some 
form  of  gain  variation  for  compensation.  For  propeller  and  early  jet 
aircraft  the  gain  variation  was  accomplished  by  a  schedule  which  was  a 
function  of  some  air  data  parameter  such  as  dynamic  pressure,  altitude, 
or  Mach  number.  In  space  vehicles  and  high-performance  supersonic  air¬ 
craft  the  gain  scheduling  becomes  very  complex  or  impossible  for  the 
following  reasons: 

1.  The  dynamic  and  static  characteristics  of  the  airframe  change 
rapidly  over  wide  ranges, 

2.  Air  data  measurements  are  not  accurate  or  are  not  available  at 
high  speeds  and  altitudes. 

3.  The  airframe  dynamic  characteristics  are  not  known  with  suf¬ 
ficient  certainty  to  permit  scheduling. 

To  eliminate  the  requirement  for  scheduling  and  air  data  measurements, 
automatic  control  systems  are  being  developed  which  include  a  self  adap¬ 
tive  loop.  The  purpose  of  the  self  adaptive  loop  is  to  maintain  optimum 
response  of  the  airframe  by  compensating  for  the  changes  in  the  dynamic 
characteristics  without  the  necessity  for  external  measurements. 
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The  Dither  Self-Adaptive  System  is  one  method  for  obtaining  optimum 
response  by  self  compensation,  A  study  of  this  system  for  pitch  control 
was  undertaken  with  the  following  objectives: 

1.  To  determine  the  ability  of  the  system  to  correct  for  the 
varying  dynamic  characteristics  of  a  flexible  airframe, 

2.  To  determine  the  optimum  variation  of  the  damping  ratio  and 
natural  frequency  of  the  short  period  mode  during  a  typical 
flight  and  how  this  variation  is  affected  by  airframe  bending 
and  by  changes  in  the  center  of  gravity  of  the  airframe. 

The  first  objective  was  investigated  by  the  comparison  of  digital  and 
analog  computer  results.  The  second  objective  was  investigated  primar¬ 
ily  from  an  empirical  approach  although  some  analytical  verification  is 
included. 
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CHAPTER  2 


DITHER  SELF  ADAPTIVE  CONTROL  SYSTEM 

2 , 1  Self  Adaptive  Systems 

Any  closed  loop  control  system  has  in  effect  some  self-adaptive 
ability*  In  the  usual  control  system  a  desirable  static  output  is  obtain¬ 
ed  under  varying  input  and  open  loop  parameter  conditions  by  proper  design 
of  the  feedback  loops.  The  name  nsel f-adapt iveH  (l)s  however*  is  reserved 
for  only  those  control  systems  which  maintain  optimum  dynamic  response 
under  varying  conditions  instead  of  the  static  response  that  normal  feed¬ 
back  optimalizes. 

Self-adaptive  systems  can  be  classified  according  to  the  ability  to 
automatically  compensate  for  either  changes  in  the  system  input  or  changes 
in  the  system  parameters*  such  as  environmental  variations.  The  ideal  self- 
adaptive  system  would  compensate  for  both  input  and  system  parameter  changes. 
Flight  control  systems  are  generally  concerned  with  parameter  changes  and 
hence,  are  designed  to  compensate  for  only  those  changes.  Since  this  in¬ 
vestigation  is  concerned  with  a  flight  control  system*  only  the  ability  of  the 
system  to  optimalize  the  dynamic  response  under  varying  aerodynamic  para¬ 
meters  will  be  considered. 

There  are  three  basic  operations  which  an  adaptive  loop  must  perform: 

1)  A  continuous  measurement  of  system  dynamic  performance  -  The 
dynamic  performance  measuring  method  should  have  a  negligible 
effect  on  the  system  response  signal  and*  in  like  manner*  the 
system  response  signal  should  have  negligible  effect  on  the 
measuring  method. 
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2)  A  continuous  evaluation  of  the  dynamic  performance  on  the 
basis  of  some  predetermined  criterion  -  The  selection  of  the 
evaluation  criterion  depends  on  the  physical  capabilities  of 
the  system.  Factors*  such  as  realizability,  complexity,  and 
cost,  may  prevent  incorporating  into  the  system  the  ability  to 
maintain  a  constant  dynamic  response.  Compromises  would  have 
to  be  made  so  that  the  system  would  be  physically  capable  of 
maintaining  optimum  a  characteristic  of  the  response;  as,  the 
relative  stability  -  damping  ratio  -  or  the  natural  frequency. 

A  figure  of  merit  -  which  is  a  number  or  a  method  of  compari¬ 
son  -  is  then  selected  on  the  basis  of  accurately  representing 
the  selected  dynamic  response  characteristic  to  be  maintained 
constant . 

3)  A  continuous  readjustment,  based  on  the  measured  performance 
and  evaluation,  of  system  control  parameters  for  optimum 
operation  --  The  methods  used  to  adjust  the  parameters  to 
obtain  optimum  operation  can  be  listed  under  three  headings: 
High  Gain  Linear  Feedback,  Programmed  Compensation,  and  Compen¬ 
sation  Using  an  Optimalizing  Controller. 

In  the  High  Gain  Linear  Feedback  Method  of  adjusting  performance  under 
varying  parameters  the  forward  loop  gain  is  maintained  high  by  compensa¬ 
tion  or  with  a  relay.  The  desired  dynamic  resporse  characteristics  are 
contained  in  the  feedback  loop.  The  high  gain  with  feedback  makes  the 
varying  parameters  in  the  forward  loop  negligible  while  the  feedback  loop 
determines  the  dynamic  response  of  the  system.  This  is  not  a  new  method  and 
is  used  in  many  electronic  applications  for  stabilization. 
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he  Pr. gtdi.Tied  ion  Method  measures  the  varying  parameter, 

or  the  condi t  ions  vhi*>  vary  the  parameter*  and  adjusts  the  compensating 
parameters  in  a  controller  in  the  forward  loop  to  counteract  the  system 
variation.  The  adjustment  is  accomplished  through  a  program  which  re¬ 
quires  a  knowledge  of  the  relationship  between  the  system  and  compensa¬ 
tion  parameters. 

The  Compensation  Using  an  Optimalizing  Controller  Method  is  a  type 
of  feedback  control  which  adjusts  parameters  so  that  a  characteristic  of 
the  system  response  Is  optimalized;  such  as5  the  damping  ratio.  In  the 
Controller  Method  an  error  signal  using  the  output  response  level  of  the 
s /stem  can  not  be  used  as  an  error  signal  for  the  parameter  adjustment  as  in 
a  normal  feedback  method  because  the  Controller  is  optimalizing  a  character¬ 
istic  of  the  output  and  not  necessarily  the  output  level.  The  system  may 
be  used  in  conjunction  with  a  model  which  gives  the  optimum  characterist ic 
in  its  output  for  a  particular  system  input.  The  Controller  samples  the 
characteristic  in  the  system  output  and  compares  the  value  with  the  optimum 
output  of  the  model.  The  Controller  then  adjusts  parameters  until  the  model 
and  system  characteristic  coincide, 

2.2  Description  of  Pitch  Control  System  to  he  Investigated. 

The  Dither  Self-Adaptive  System  is  capable  of  controlling  a  vehicle 
about  its  pitch*  roll*  and  yaw  axes.  This  investigation*  however*  will 
be  concerned  only  with  pitch  control.  The  system  which  is  to  be  investi¬ 
gated  is  shown  in  Figure  2,1.  The  position,  rate  and  acceleration  sensors 
that  are  required  for  the  feedback  loops  are  assumed  to  have  a  unity  trans¬ 
fer  function.  The  rectifier  in  the  self  adaptive  loop  is  also  assumed  to 
have  a  unity  transfer  function. 
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INVERSE  MODEL 


Figure  2.1.  Dither  System 


In  the  conventional  loop  an  error  signal,  S  c,  is  formed  by  comparing 
the  feedback  signal  to  the  command  input  signal*  The  error  signal  is 
amplified  by  the  adaptive  gain,  kg-  ,  and  the  resulting  signal  operates  the 
servo  actuator*  This  investigation  involves  a  missile  which  has  gimbaled 
thrust  chambers*  The  servo  actuator  controls  the  direction  of  the  chamber 
axis  which  results  in  a  thrust  vector  signal  <5*t.  The  thrust  vector  signal 
affects  the  vehicle  performance  which  is  represented  by  the  block  labeled 
"Vehicle  Dynamics"*  The  reaction  of  the  vehicle  airframe  to  the  change  in 
direction  of  the  thrust  axis  is  detected  by  position,  rate,  and  acceleration 
sensors  and  is  fed  back  through  the  Inverse  Model  to  the  input*  The  function 
of  the  Inverse  Model  will  be  explained  in  Section  2.3. 

For  the  self  adaptive  loop  a  dither  signal  is  introduced  at  the  input* 

This  signal  is  then  amplified  by  kg-  ,  operates  the  servo  actuator,  and 

* » 

excites  the  vehicle  dynamics.  The  pitch  acceleration,  0,  is  passed  through 
a  bandpass  filter  where  the  dither  frequency  is  separated  from  the  command 
input  signal  frequencies.  The  dither  frequency  is  then  rectified  and  com¬ 
pared  with  a  reference  signal.  The  error,  ;  operates  an  integrating 
servo  which  varies  the  adaptive  gain,  k g  ,  until  the  reference  and  rectified 
signals  are  equal. 

2 . 3  Principles  Used  to  Maintain  Optimum  Dynamic  Response. 

When  a  physical  aerodynamic  parameter  changes,  all  the  dynamic 
characteristic  parameters  of  the  system  are  affected.  Typical  physical 
aerodynamic  parameters  are  Mach  number  and  altitude,  and  typical  dynamic 
characteristic  parameters  are  natural  frequencies  and  damping  ratios.  In 
terms  of  the  system  transfer  function  a  change  in  aphysical  parameter  could 
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shift  the  posit  ion  of  all  the  poles  and  zeros  of  the.  transfer  function* 

If  constant  dynamic  response  was  desired  for  the  system,  the  variation  in 
all  the  poles  and  zero  would  have  to  be  counter-balanced*  The  complexity 
of  the  situation  can  be  visualized  by  placing  a  root  locus  plot  alongside 
a  missile  traveling  many  times  the  speed  of  sound*  The  physical  parameters 
of  the  missile  are  constantly  changing,  and  each  parameter  would  be  varying 
all  the  poles  and  zeros  on  the  root  locus.  A  typical  root  locus  with  vary* 
ing  poles  and  zeros  is  shown  in  Figure  2*2.  A  controller  would  have  to  be 
built  to  sense  each  of  the  variations  and  then  decide  how  to  counter* 
balance  by  compensation.  In  practice  exact  compensation  is  too  complicated 
to  be  useful  and  approximations  are  always  made. 

Simplicity  is  one  of  the  main  objectives  of  the  dither  adaptive  loop. 

The  only  components  that  have  been  added  to  the  existing  ptich  control 
system  are  a  dither  generator,  filter,  rectifier,  and  a  comparison- integra¬ 
tion  network.  No  variable  shaping  networks  are  used,  and  all  adjustable 
compensation  is  accomplished  by  varying  the  gain  of  the  system. 

Aviators  are  primarily  affected  physically  by  the  short  period  mode 
of  the  pitch  aerodynamic  equations.  Tests  have  indicated  that  pilots  prefer 
a  damping  ratio  of  0.7  and  a  natural  frequency  of  3  radians  per  second  for 
this  mode  (2).  Since  the  missile  used  in  this  investigation  is  designed 
for  manned  flights,  the  damping  ratio  of  0.7  and  frequency  of  3  were  de¬ 
sired.  Since  the  dither  system  does  not  use  variable  shaping  networks,  it  is 
physically  impossible  to  maintain  both  the  damping  ratio  and  frequency  con¬ 
stant  at  the  desired  values.  Therefore,  a  damping  ratio  of  0.7  with  a 
minimum  natural  frequency  variation  are  defined  as  the  optimum  dynamic 
characteristics  of  the  short  period  mode  during  a  flight. 
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There  are  actually  two  methods  being  employed  In  the  overall  system 
for  varying  the  control  parameter  to  obtain  optimum  performance*  The 
control  parameter  of  the  system  is  the  gain*  The  two  methods*  which  are 
discussed  in  Section  2.1,  are  High  Gain  Linear  Feedback  and  Programmed 
Compensat ion. 

If  the  gain  of  the  system  could  be  maintained  at  a  high  value,  the 

desired  dynamic  response  characteristic  could  be  inserted  in  a  feedback 

loop  and  would  determine  the  dynamic  response  of  the  overall  system. 

This  is  the  principle  of  the  High  Gain  Linear  Feedback  Method  shown  in 

Figure  2.3.  H(s)  is  referred  to  as  the  "Inverse  Model'1  since  the  inverse 

of  its  transfer  function  determines  the  dynamic  response  of  the  system. 

If  a  damping  ratio  of  0.7  and  a  natural  frequency  of  3  radians  per  second 

2 

were  desired,  the  Inverse  Model  would  have  a  transfer  function  of  s  +4.2s+9. 
Unfortunately,  the  varying  dynamic  characteristics  of  the  airframe,  as 
represented  by  G2(s),  cannot  be  submerged  by  merely  using  high  gain  feed* 
back  loops  (3).  The  required  aerodynamic  control  moments  and  power  required 
for  high  gains  in  flight  control  systems  are  not  available  due  to  limits 
imposed  by  vehicle  design  considerations  unrelated  to  the  automatic  control 
system  design.  Even  if  high  gains  could  be  obtained,  the  gains  could  re¬ 
sult  in  severe  bending  and  slosh  mode  oscillations  when  considering  flexible 
airframes.  The  values  that  the  Inverse  Model  transfer  function  can  assume 
and  still  affect  the  short  period  mode  in  the  desired  manner  are  also  limited 
for  flexible  airframes.  The  limit  on  the  values  of  the  transfer  function 
for  this  investigation  is  discussed  in  Section  2.5. 

For  flight  control  systems  the  limitation  on  high  gains  and  on  values 
of  the  Inverse  Model  prevents  the  use  of  a  true  High  Gain  Linear  Feedback 
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G2(s)'  =  >^G2*(s) 


If  is  .large 
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where  is  the  aerodynamic  gain. 

K,fGi  (s)Gp(s)  _ 
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and  K^H(s)Gx(sJG2(s)  >}  1, then 
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©c  H(s). 


figure  2.3.  General  feedback  System. 


self  ad-ptive  scheme 


Inverse  Model  can  be  used  as  a  fixed 


'  1  \  ir  >  f  tf, 

shaping  network  to  improve  the  dynamic  characteristics  of  the  short  period 
mode  and  is  employed  as  such  in  the  system  under  investigation* 

In  view  of  restricting  the  Inverse  Model  to  the  idea  of  a  fixed  shaping 
network,  the  self  adaptive  scheme  used  in  the  Dither  System  can  be  thought 
of  as  a  form  of  the  Programmed  Compensation  Method*  The  three  basic  opera¬ 
tions  that  an  adaptive  loop  must  perform  are  listed  in  Section  2*1.  The 
manner  in  which  the  Dither  System  performs  these  operations  are  listed  below: 

1,  The  dither  signal  is  the  means  for  the  continuous  measurement 
of  the  system  dynamic  response*  In  order  to  have  a  negligible 
effect  on  the  system  response  signal,  the  dither  frequency 
should  be  five  or  more  times  the  short  period  mode  frequency 
and  the  input  level  should  be  small  so  that  the  oscillations 
are  not  noticeable  to  the  vehicle  personnel.  The  reason  for 
using  pitch  acceleration  instead  of  pitch  angle  or  rate  as  the 
pickoff  for  the  adaptive  loop  is  to  enable  smaller  inputs  for 
the  same  output.  The  upper  limit  for  the  dither  frequency 
depends  on  the  ability  of  the  actuator  to  respond  to  the  fre¬ 
quency  and  on  the  bandwidths  of  the  rest  of  the  system  components. 
The  upper  limit  is  around  40  radians  per  second  for  available 
systems  . 

2.  The  continuous  evaluation  of  the  dynamic  performance  is  per¬ 
formed  by  comparing  the  amplitude  of  the  pitch  acceleration  to 
dither  input  level  evaluated  at  the  dither  frequency  to  the 
Performance  Criterion  Reference.  It  remains  to  be  discussed 
below  and  in  Chapter  4  how  accurately  the  dynamic  performance 
is  being  evaluated.  In  this  investigation  the  damping  ratio 
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is  being  used  as  the  measure  of  dynamic  performance, 

3.  The  continuous  readjustment  of  the  system  control  parameter, 
which  is  the  gain,  is  accomplished  by  a  servo  driven  by  the 
error  resulting  from  the  evaluation  procedure* 

In  the  Programmed  Compensation  Method  the  parameter  adjustment  is  ac¬ 
complished  through  a  program  which  requires  a  knowledge  of  the  relationship 
between  the  system  and  compensation  parameters.  In  the  dither  system  the 
program  consists  of  a  single  reference  level  which  requires  a  knowledge 
of  the  variation  of  the  system  damping  ration  with  the  amplitude  of  the 
acceleration  to  dither  input  level  ratio* 

The  degree  of  self  adaptability  of  the  dither  system  depends  on  the 
accuracy  with  which  the  damping  ratio  is  measured  by  the  acceleration  to 
input  level  ratio*  There  are  many  figures  of  merit  that  can  be  used  as  a 
criterion  for  the  dynamic  performance  of  a  system*  Some  of  the  figures  of 
merit  are 

is  an  error  signal*  The  applicability  of  these  figures  of  merit  for  use 
in  establishing  a  damping  ratio  of  0.7  has  been  determined  for  a  second 
order  system  to  be  as  follows:  (4) 

oO 


^6  selects  zero  damping  ratio  as  optimum* 

selects  0.5  damping  ratio  as  optimum* 

S  l£t  selects  0.7  damping  ratio  as  optimum. 

o 

selects  zero  damping  ratio  as  optimum. 

^ selects  0*7  damping  ratio  as  optimum. 

If  the  system  was  of  second  order,  both  /£/  and  ^ 

would  select  0.7  as  the  optimum  damping  ratio  based  on  the  figure  of 
merit  approaching  a  minimum  value  at  a  damping  ratio  of  0.7.  The 

figure  of  merit  has  the  sharpest  minimum  and  would  be  the 
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preferred  criterion.  With  a  flexible  airframe  the  system  can  not  usually 
be  approximated  by  a  second  order  system  with  any  accuracy,  and  the  ap¬ 
plicability  of  the  above  criterion  for  establishing  optimum  damping  ratics 
of  the  short  period  mode  under  these  conditions  is  not  known.  The  corn- 


criterion  into  a 


ponents  required  to  mechanize  the 


system  is  certainly  more  complex  than  the  Dither  System,  and  the  possibility 
of  the  criterion  selecting  a  damping  ratio  of  0.7  for  a  non  second  order 
system  is  considered  doubtful.  A  different  figure  of  merit  called  the 
Impluse  Response  Area  Ration  (IRAR)  will  permit  the  selection  of  any  de¬ 
sired  damping  ratio  and  have  the  error  signal  assume  a  zero  value  when  the 
response  is  at  the  selected  damping  ratio  (5).  The  method  required  for  use 
of  the  IFAR  criterion  would  be  to 'pulse  the  system  and  then  measure  the 
area  of  the  response.  The  area  is  a  measure  of  the  damping  ratio.  By 
comparing  the  area  with  a  reference  area  corresponding  to  the  desired 
damping  ratio  an  error  signal  is  developed  which  would  vary  the  gain  of 
the  system  until  the  error  was  zero.  In  a  flexible  airframe  control  system 
which  does  not  use  shaping  networks,  the  poles  and  zeros  of  the  vehicle 
transfer  function  are  continually  varying.  For  the  reference  area  to  in¬ 
dicate  the  true  damping  ratio  it  would  also  have  to  vary  with  the  parameters 
and  in  effect  would  need  a  self  adaptive  loop  to  adjust  its  value.  The 
system  would  be  very  complex  if  not  impossible  to  mechanize. 

The  Dither  System  works  on  the  principle  of  maintaining  the  gain  of  the 
system  constant  throughout  the  flight  and  assumes  that  a  constant  gain  will 
maintain  an  optimum  damping  ratio.  Referring  to  Figure  2.3  and  substituting 


Af/p  25  Qf{  an  adaptive  loop  is  formed. 


Or  (M  u/0)  z  fa  Ms  (kg  — . 

/+  /%  Ms  G,  (j'UA)  6\»( ) 
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Since  Ks  i%  &i  &&  )  has  a  denominator 

/“>  tl 

of  greater  order  than  the  numerator ,  ^  Hi  M<£  &/  (^^c)  when 

Wo  is  a  large  value.  Therefore,  is  a  direct  measure  of 

the  gain,  »  of  the  system.  The  self  adaptive  loop  maintains  the 

value  of  Kj  equal  to  the  reference  value.  If  a  constant  value  of 

K^Ms  an  accurate  measure  of  the  damping  ratio,  £  ,  for  a  flight 

control  system,  then  a  plot  of  R§,M$  versus  £  should  show  a  constant 
RcpA&  for  a  particular  £  for  all  flight  cases.  Using  the  R\s  !f  of  the 
flight  cases  to  be  investigated  a  plot  of  KfMg  vs  ^  is  shown  in 
Figure  2.4  for  aerodynamic  equations  with  no  bending  modes  included,  the 
first  bending  mode  included  and  the  first  and  third  bending  modes  included. 
The  aerodynamic  equations  and  conditions  are  discussed  in  Section  2.4.  As 
is  observed  from  the  Figure,  ifoMf  is  an  accurate  criterion  for  a  con¬ 
stant  £  for  no  bending,  but  as  bending  modes  are  included  and  the  order 
of  the  system  increases  KgMf  becomes  a  less  accurate  figure  of  merit. 

The  investigation  of  the  variation  of  the  damping  ratio  for  the  various 
aerodynamic  conditions  and  the  determination  of  the  optimum  Performance 
Criterion  Reference  and  the  optimum  dither  frequency  is  contained  in 
Chapter  4. 
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2.4  Aerodynamic  Equations  and  Modifications* 

The  airframe  used  in  this  investigation  is  the  two  staged  Saturn 
Booster*  Perturbated  aerodynamic  equations  with  coefficients  for  four 
flight  cases  are  included  in  Appendix  A.  The  flight  cases  are  identified 
by  the  time  in  seconds  after  launch.  The  flight  cases  are  t  =  20*  t  =  75, 
t  =  101,9,  and  t  =  200,  The  t  =  200  case  is  for  the  second  stage  only, 
as  the  first  stage  separates  before  t  =  200  seconds.  The  reference  for 
the  equations  is  the  nose  of  the  second  stage.  For  flight  cases  t  =  20,  t  = 
75,  and  t  =  101.9  the  acceleration,  rate,  and  position  sensors  are  located 
on  the  same  platform  in  the  first  stage.  The  variation  in  performance  for 
various  platform  positions  will  not  be  covered  in  this  investigation.  The 
platform  position  used  had  been  selected  for  maximum  system  stability 
immediately  after  launch  and  for  airframe  considerations.  The  platform 
position  for  case  t  =  200  is  in  the  nose  of  the  second  stage.  The  transfer 
functions  which  were  determined  from  the  equations  using  digital  computers 
are  included  in  Appendix  A. 

The  aerodynamic  equations  contain  three  bending  modes  and  one  fuel 
slosh  mode.  A  typical  root  locus  is  shown  in  Figure  2.5.  It  will  be 
observed  from  the  Figure  that  the  open  loop  system  is  unstable.  Case  t  = 

20,  which  is  the  first  case  considered  after  launch  and  where  the  missile 
is  still  verticle,  is  the  only  flight  case  with  a  stable  open  loop.  An 
unstable  aerodynamic  condition  occurs  when  the  center  of  pressure  is  for¬ 
ward  of  the  center  of  gravity.  A  missile  of  conventional  shape  will  be 
unstable  unless  stabilizing  fins  are  attached  because  the  center  of 
pressure  of  the  fuselage  alone  is  close  to  the  shoulder  whereas  the  center 
of  gravity  is  close  to  the  geometric  center  (6).  For  large  missiles,  such 
as  the  Saturn,  the  trend  is  towards  stabilization  through  a  control  system 
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rather  than  providing  for  a  stable  airframe  because  of  the  weight  required 
for  the  stabilizing  fins.  The  control  torques  required  for  the  control 
system  are  provided  for  in  the  Saturn  by  gimbaled  thrust  chambers.  It 
will  also  be  observed  from  the  Figure  that  the  Slosh  Mode  has  a  negligible 
effect  on  the  system  stability.  Multitank  design  with  slosh  suppressors 
are  used  in  the  Saturn  to  obtain  this  condition.  The  different  slosh  mode 
frequencies  correspond  to  different  tanks. 

The  second  bending  mode  becomes  unstable  for  an  adaptive  variable  gain 
of  approximately  0.1  for  the  t  =  20,  75,  and  101.9  cases.  In  order  to 
eliminate  the  possibility  of  obscuring  the  actual  dither  system  performance 
by  the  complexity  of  adding  bending  mode  cancellation  schemes,  the  unstable 
second  bending  mode  was  handled  in  one  of  two  ways: 

1.  The  second  bending  mode  was  eliminated  from  the  aerodynamic 
equations  assuming  that  a  bending  cancellation  method  was  in  use. 

2.  The  second  bending  mode  was  eliminated  by  assuming  that  a 
compensation  method  was  in  use,  such  as  a  frequency  tracker  and  notch 
filter  scheme,  which  placed  a  pole  and  zero  over  the  open  loop  zero  and 
pole  respectively  of  the  bending  mode. 

A  bending  cancellation  method  which  could  be  used  for  procedure  1 


above  is  shown  in  Figure  2.6.  This  is  one  of  many  cancellation  methods 
proposed  by  the  Autonetics  Division  of  North  American  Aviation.  The  purpose 
of  the  system  is  to  form  a  signal  of  equal  magnitude  and  opposite  phase  to 
the  bending  mode  signal.  The  signal  formed  is  then  added  to  the  control 
signal  and  effectively  cancels  the  bending  mode  from  the  system  response 
signal.  Bandpass  filters  land  2  are  variable  tuned  filters  with  center 
frequencies  of  lUo  .  If  filter  1  is  not  tuned  to  the  bending  frequency, 
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Figure  2.6.  Bending  Cancellation  Method 


the  bending  frequency  i^n  1  component  is  phase  shifted  by  an  amount  (p 
and  the  output  signal  from  filter  1  is  sin  (p')  •  The  bending  sig¬ 

nal  is  again  phase  shilled  by  an  amount  (p  by  filter  2  with  output  signal 
being  sin^2 The  outputs  from  the  two  filters  pass  through 
phase  shifters  which  produce  a  90  degree  phase  shift  between  the  signals. 
The  shifted  signals  are  then  multiplied,  and  the  resulting  signal  drives 
an  integrator  which  varies  the  center  frequencies  of  the  filters  until  the 
output  from  the  multiplier  is  zero.  As  is  shown  below,  the  output  of  the 
multiplier  is  zero  when  the  filters  are  tuned  to  the  bending  frequency. 
Output  of  multiplier  =  sin  (ciu'oJ&t  cos  ( U/osZ^+  (p) 

=  (sin  cos  p  -i  cos  sxnp)  (cos  lt/os£*cos£p -sin  *U/0  t  sin  Ip). 

The  integrator  is  not  sensitive  to  frequencies  much  greater  than  the 
bending  frequency. 

Therefore, 

Effective  Output  =  1/2  cos  2  <p)  sin^  -  1/2  cos  p?  sin  2^2? 

-  1/2  sxnp 
=  0  when  p  is  zero. 

With  the  filters  tuned  to  the  bending  frequency,  the  180  degree  phase  shift 
is  accomplished  by  subtraction  from  the  feedback  signal.  The  amplitude  of 
the  bending  signal  is  adjusted  by  the  rect i fier- integrator  comparison  loop 
using  the  outputs  of  filters  l  and  2,  for  the  comparison.  Another  bend¬ 
ing  mode  cancellation  method  using  a  rejection  filter  philosophy  is  con¬ 
tained  in  Reference  7. 

The  different  aerodynamic  conditions  that  are  investigated  are  listed 
below: 

1.  No  bending  and  no  slosh  modes  included  In  the  equations. 

The  resulting  equations  are  the  rigid  body  equations. 
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2.  The  fits*  bending  mode  added  to  the  rigid  body  equations* 

3*  The  first  and  third  bending  modes  added  to  the  rigid  body 
equations*  The  second  bending  mode  is  assumed  removed  by  a  bending  mode 
cancellation  method* 

4.  The  first  and  third  bending  and  first  slosh  modes  added  to 
the  rigid  body  equations*  The  second  bending  mode  is  assumed  removed  by 
a  frequency  tracker  and  notch  filter  method* 

In  Chapter  3  the  analog  computer  investigation  was  conducted  using 
aerodynamic  conditions  ls  2*  and  3  and  flight  cases  t  =  20*  75*  101*9  and 
200.  The  second  bending  mode  was  included  in  condition  3  for  flight  case 
t  =  200  because  the  mode  is  not  unstable  for  this  case* 

In  Chapter  4,  the  empirical  investigation  was  conducted  using  aero¬ 
dynamic  conditions  ls2,3,  and  4  and  flight  cases  t  =  20*  75*  and  101*9* 

2 . 5  Limitations  on  Inverse  Model  by  Flexible  Airframe* 

The  system  under  investigation  is  unstable  unless  a  shaping  network 
is  employed.  From  a  root  locus  viewpoint  if  a  complex  pair  of  zeros  were 
provided  in  the  left  half  plane*  the  unstable  short  period  mode  would  be 
shaped  into  the  left  half  plane  and  the  system  would  be  stabilized*  The 
Inverse  Model  provides  these  shaping  zeros* 

If  a  damping  ratio  of  0.7  and  a  natural  frequency  of  3  radian  per 
second  were  the  desired  dynamic  characteristics  of  the  short  period  mode* 
the  Inverse  Model  should  ideally  provide  a  pair  of  zeros  at  s  -  -2*1  _  j 2 . 18 
on  the  root  locus.  A  high  system  gain  would  then  establish  a  complex  pair 
of  roots  of  the  system  characteristic  equation  at  values  s  -  -2.1+j  2*18. 

The  roots  would  fix  the  dynamic  characteristics  of  the  short  period  mode 
at  a  damping  ratio  of  0*7  and  a  frequency  of  3*  As  was  discussed  in  Section 
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2.3,  high  g^m  capabilities  an'  not  available  in  flight  control  systems. 

The  desired  dynamo  characteristics  could  still  be.  obtained  even 
though  high  gains  are  not  available  by  increasing  the  real  and  imaginary 
coordinates  of  the  Inverse  Model  zeros.  The  zeros  could  shape  the  short 
period  dynamic  characteristic  curve  in  such  a  manner  that  for  lower  values 
of  gain  the  desired  characteristics  are  available.  Increasing  the  coordin¬ 
ates  of  the  Inverse  Model  zeros  is  possible  in  a  rigid  airframe  with 
negligible  center  of  gravity  movement.  For  flexible  airframes  the  bend¬ 
ing  modes  must  be  considered  in  establishing  the  position  of  the  Inverse 
Model  zeros. 

If  the  coordinates  of  the  zeros  are  increased  beyond  certain  values 
with  a  flexible  airframe  system,  the  zeros  thape  primarily  a  different 
mode  than  the  short  period  mode.  While  the  short  period  mode  would  be 
stabilized,  its  maximum  damping  ratio  and  natural  frequency  values  would 
be  less  and  its  coupling  with  the  shaped  mode  would  be  gi eater  than  when 
the  zeros  shape  primarily  the  short  period  mode.  A  digital  computer  investi¬ 
gation  of  the  optimum  position  of  the  zeros  was  conducted. 

The  preliminary  investigation  showed  that  flight  case  t  -  75  was  the 
critical  flight  case  and  that  the  two  bending  modes  plus  slosh  mode  condi¬ 
tion  was  the  critical  condition.  Flight  case  t  -  75  was  further  investi¬ 
gated  as  to  the  benefit  of  using  rate  and  position  feedback  which  would 
give  one  real  axis  zero  or  of  using  acceleration,  rate9  and  position  feed¬ 
back  which  would  give  the  complex  pair  of  zeros.  The  complex  zero  arrange¬ 
ment  provided  for  more  shaping  of  the  short  period  mode  and  less  possible 
coupling  between  the  first  bending  mode  and  short  period  mode.  The  coordin¬ 
ates  of  the  complex  zeros  were  then  varied  to  obtain  optimum  values. 
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Typical  root  loci  are  shown  in  Figure  2.7.  For  Inverse  Model  position  3, 
the  zeros  shape  primarily  the  first  bending  mode.  While  the  short  period 
mode  is  stabilized,  it  is  not  being  shaped  in  an  optimum  manner  and  has 
the  disadvantages  of  decreasing  in  damping  ratio  as  gain  increases  beyond 
a  certain  value  and  of  being  coupled  with  the  bending  mode.  The  short  period 
mode  could  never  obtain  a  0.7  damping  ratio.  For  Inverse  Model  positions 
1  and  2,  the  short  period  mode  is  being  primarily  shaped.  The  bending  mode 
is  also  shaped  further  into  the  left  half  plane  increasing  its  stability 
while  producing  a  negligible  increase  in  coupling.  Position  2,  with  the 
zeros  at  s  =  1.8  +j. 6,  appears  to  be  the  optimum  position.  The  maximum 
natural  frequency  of  the  short  period  mode  at  a  damping  ratio  of  0.7  occurs 
at  this  position.  As  is  observed,  the  natural  frequency  will  never  reach 
3  radians  per  second.  A  more  complicated  shaping  network  could  be  used  to 
increase  the  natural  frequency,  but  since  one  of  the  principal  advantages 
of  the  Dither  System  is  simplicity,  a  compromise  has  been  made  on  the  desired 
dynamic  characteristics.  As  was  defined  in  Section  2.3,  the  optimum  dynamic 
characteristics  of  the  short  period  mode  for  the  investigation  are  a  damp¬ 
ing  ratio  of  0.7  with  a  minimum  natural  frequency  variation. 

The  optimum  Inverse  Model  transfer  function  for  the  severest  flight 

2 

case  and  condition  will  remain  fixed  at  s  -»  3. 6s+3.6  for  all  flight  cases 
and  for  all  aerodynamic  conditions. 
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Figure  2.7.  Variation  of  Inverse  Model 
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CHAPTER  3 


COMPUTER  SIMULATION 


3  .  1  Int.  rod  jcLi.uTi . 

The  object  of  this  Chapter  is  to  investigate  the  ability  of  the 
Dither  Self  Adaptive  System  to  correct  for  the  varying  dynamic  character- 
lstics  of  a  flexible  airframe  during  a  typical  flighto  The  correction  is 
made  by  the  adjustment  of  the  control  parameter  which  is  the  adaptive  gain, 

. 

An  analog  computer  simulation  was  conducted  for  the  primary  purpose 
of  determining  the  values  of  the  adaptive  gain  which  the  system  selected 
for  the  various  flight  cases  and  for  the  various  aerodynamic  conditions. 

The  values  of  gain  which  were  selected  in  the  computer  simulation  are 
compared  with  the  values  of  gain  determined  by  digital  computer  methods. 
Digital  computer  techniques  are  employed  in  the  investigation  in  Chapter  4, 
and  the  comparison  with  analog  computer  results  in  this  Chapter  is  used  as 
an  indication  of  the  accuracy  of  the  results  of  Chapter  4  when  the  transi- 
tion  is  made  from  numerical  to  electrical  simulation.  The  damping  ratios 
of  the  short  period  modes  which  are  associated  with  the  gains  selected  by 
the  analog  simulation  for  the  various  aerodynamic  conditions  are  determined 
from  root  loci  using  digital  computer  techniques  and  are  compared  to  show 
the  effect  of  a  flexible  airframe  on  the  ability  of  the  dither  system  to 
maintain  a  constant  damping  ratio.  The  damping  ratios  of  the  short  period 
mode  resulting  from  the  use  of  the  dither  system  are  also  compared  with  the 
damping  ratios  resulting  from  the  use  of  a  constant  K£  system. 

The  dither  frequency  and  reference  level  which  would  maintain  through 
the  flight  an  average  damping  ratio  of  the  short  period  mode  of  0.7  with  a 
minimum  variation  of  the  damping  ratio  and  natural  frequency  were  not  used 
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in  the  c.  input  cr  situ! -it  i»n»  The  icsults  of  the  computer  simulation 
indicated  tbu  necessity  ri  determining  optimum  frequencies  and  reference 
levels.  The  investigation  of  this  area  is  contained  in  Chapter  4. 

The  secondary  purpose  of  the  analog  computer  simulation  was  to  deter¬ 
mine  for  the  various  flight  cases  and  aerodynamic  conditions  the  rate  at 
which  the  adaptive  gain  was  adjusted  and  the  reaction  of  the  self  adaptive 
loop  to  wind  gusts,  actuator  movements*  and  loss  of  the  dither  signal.  A 
complete  analysis  of  the  self  adaptive  loop  was  not  conducted  in  the  investi¬ 
gation;  and  therefore,  the  values  of  the  loop  gain  and  filter  bandwidth  for 
optimum  dynamic  response  of  the  loop  were  not  determined.  The  investigation 
was  conducted  using  different  loop  gains  for  the  various  flight  cases.  A 
qualitative  discussion  of  the  adaptive  loop  is  contained  in  Section  3.5. 

3. 2  Analog  Computer  Setup. 

The  analog  investigation  was  conducted  using  aerodynamic  conditions 
1,2,  and  3  as  described  in  Section  2.4.  For  Flight  Case  t  =  200  the 
second  bending  mode  was  included  because  the  mode  is  not  unstable  for 
this  case.  The  system  was  simulated  in  two,  forty  amplifier.  Pace  Computers 
by  Electronic  Association  Inc.  The  voltmeter  used  to  determine  the  K  <$ 
values  and  the  integrating  servo  and  rectifier  used  in  the  adaptive  loop 
were  inherent  components  of  the  computers.  Two  six  channel  Sanborn  re¬ 
corders  were  used.  A  Hewlett-Packard,  Low  Frequency,  Function  Generator  was 
used  as  the  source  of  the  dither  signal. 

The  computer  simulation  diagram  is  shown  in  Figure  B-l.  The  reference 

%  •  * 

values  of  pitch  acceleration,  0,  pitch  rate  0,  and  pitch  angle,  0,  which 
are  measured  at  the  nose  of  the  second  stage  are  converted  to  the  platform 

4  •  » 

values,  9p  ,  0p  ,  and  Gp  ,  for  flight  cases  t  =  20,  75,  and  101.9,  be¬ 
cause  the  sensors  for  these  cases  are  on  the  same  platform  which  is  located 
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In  the  first  stage.  The  reference  values  are  the  platform  values  for 
Flight  Case  t  -  200,  because  the  platform  for  this  case  is  in  the  nose  of 
the  second  stage. 

The  first  flight  case  simulated  was  t  =  75  for  aerodynamic  condition 
3.  Static  and  dynamic  checks  were  conducted.  The  procedure  used  for  the 
dynamic  checks  was  to  excite  the  bending  and  short  period  modes  by  voltage 
steps  introduced  into  the  amplifiers  of  the  various  mode  simulations  and 
to  compare  the  resulting  frequencies  and  damping  ratios  with  the  frequencies 
and  damping  ratios  determined  by  digital  computer  techniques.  The  trans¬ 
fer  function  resulting  from  the  removal  of  the  poles  and  zeros  associated 
with  the  second  bending  mode  from  the  vehicle  dynamic  transfer  function, 
which  was  derived  from  the  aerodynamic  equations  with  three  bending  modes 
and  one  slosh  mode,  was  considered  a  good  approximation  for  determining  the 
first  and  third  bending  mode  frequencies  for  condition  3,  A  comparison  with 
the  analog  results  of  the  dynamic  checks  indicated  that  new  transfer  func¬ 
tions  based  on  the  applicable  aerodynamic  equations  were  required.  The 
transfer  function  for  flight  case  t  =  75  for  aerodynamic  condition  3  in¬ 
dicated  that  the  only  poles  and  zeros  of  the  transfer  function,  which  was 
derived  from  the  aerodynamic  equations  with  three  bending  modes  and  one 
slosh  mode,  that  were  appreciably  moved  were  the  zeros  associated  with  the 
third  bending  mode.  The  value  of  the  zeros  were  changed  from  s  =~11.56+j21, 
63  to  s  =  -.17+  J24.3.  The  effect  of  the  bending  and  slosh  mode  parameters 
on  the  vehicle  dynamic  characteristics  can  be  determined  by  observing  the 
movement  of  the  poles  and  zeros  of  the  transfer  functions,  which  are  contain 
ed  in  Appendix  A,  for  the  various  aerodynamic  conditions.  The  static  checks 
agreed  to  within  one  half  of  a  per  cent  with  analytical  values.  The  mode 
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frequencies  of  the  analog  computer  agreed  to  within  two  per  cent  of  the 
frequencies  determined  by  digital  methods,  and  the  damping  ratios  agreed 
to  within  five  per  cent. 

The  procedure  used  to  simulate  the  self  adaptive  loop  was  as  follows: 

1.  The  filter  was  simulated  and  the  bandwidth  was  varied  until 
a  sharp  peaked  frequency  response  was  obtained  with  the  peak  at  the  dither 
frequency  of  30  radians  per  second.  A  bandwidth  of  3  radians  per  second 
produced  a  sharp  peak. 

2.  Two  diodes  which  were  inherent  components  of  the  computer 
were  used  to  form  a  rectifier.  The  outputs  from  the  two  diodes  were  of 
different  values;  and  therefore, gain  pots  were  included  with  each  of  the 
diodes  so  that  equal  outputs  could  be  obtained. 

3.  The  adaptive  gain,  k£  ,  was  simulated  by  using  a  linear 
gain  pot  driven  by  an  integrating  servo/  The  servo  and  pot  were  inherent 
components  of  the  computer.  One  volt  into  the  servo  resulted  in  a  gain 
change  of  0.2,  A  linear  pot  was  used  because  it  was  readily  available. 

A  shaped  gain  pot  should  be  used  in  reality  to  obtain  a  small  gain  to  servo 
input  ratio  for  low  values  of  gain  and  a  large  ratio  for  large  values  of 
gain.  A  small  gain  to  servo  input  ratio  is  desired  at  low  values  of  gain 
so  that  a  wind  gust  or  other  disturbance  would  not  reduce  the  gain  apprec  - 
ably  with  the  possibility  of  the  system  going  unstable.  As  the  gain  in¬ 
creases  the  dynamics  of  the  short  period  mode  improve  so  that  while  a 
disturbance  at  high  gains  would  vary  the  gain  more  than  at  low  gains  the 
response  time  is  improved  and  stability  considerations  are  not  as  critical. 
Compromises  must  be  made  in  the  shaping  of  the  gain  pot  between  the  large 
ratio  of  gain  to  servo  input  desired  for  fast  adaptive  response,  the  small 
ratio  desired  when  disturbances  are  experienced,  and  the  ratio  required  for 
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constant  dynamic  characteristics  of  the  adaptive  loop  itself  during 
a  flight.  The  discussion  of  maintaining  the  adaptive  loop  gain  appreciably 
constant  during  a  flight  is  contained  in  Section  3.5. 

4.  A  gain  pot  was  inserted  before  the  integrating  servo  to 
provide  for  varying  the  ratio  of  to  servo  input  and  the  adaptive  loop 

gain.  The  adaptive  loop  gain  was  adjusted  by  varying  the  gain  pot*  K  A  , 

until  the  gain  adjustment  response*  K£  ,  ceased  to  oscillate  and  was 
critically  damped  with  a  suitable  rate  of  change  of  the  value  of  . 

The  Reference  DC  signal  was  then  adjusted  to  provide  for  a  value  of 
of  1.  The  adjustment  of  the  adaptive  loop  gain  to  give  a  suitable  rate  of 
change  of  and  the  setting  of  the  Reference  level  so  that  had  a  value 

of  1  was  performed  for  each  aerodynamic  condition  of  Flight  Case  t  ~  75 
in  order  to  obtain  reference  values.  The  reference  values  that  were  obtain¬ 
ed  eliminated  the  necessity  for  a  complete  evaluation  of  the  analog  simulat¬ 
ed  loop  and  were  used  in  the  other  flight  cases  to  form  the  basis  for  the 
comparison  of  the  values  of  I<£  that  the  system  selected. 

When  the  system  was  energized,  the  servo  chattered  because  it  was 
following  the  half  cycle  sinusoidal  like  waveform  at  its  input.  A  smooth¬ 
ing  network  would  be  required  at  the  rectifier  output  to  produce  a  flatter 
output  and  eliminate  the  servo  chatter.  In  order  to  maintain  the  adaptive 
loop  as  simple  as  possible  and  also  to  investigate  the  effect  of  bending 
interference  the  bandwidth  of  the  filter  was  increased  to  5  radians.  The 
third  bending  mode  frequency  which  was  now  at  the  limits  of  the  bandwidth 
was  passing  through  the  loop  because  the  waveform  from  the  rectifier  flat¬ 
tened  out  so  that  the  servo  stopped  chattering.  Since  a  dither  frequency 
of  30  radians  per  second  excited  the  third  bending  mode  for  aerodynamic 
condition  3,  the  computer  results  would  indicate  the  performance  of  the 
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ten  when  bending  or  noise  signals  were  contained  in  the  filter¬ 
ed  ft  ti  cy  of  the  self  adaptive  loop.  In  order  to  determine  the  per- 
formarc  ^  i  the  system  with  no  bending  and  minimum  noise  signals  in  the 
self  adaptive  loop,  reference  values  were  also  recorded  for  Flight  Case 
t  =  75  for  aerodynamic  condition  3  with  a  dither  frequency  of  20  radians 
per  second.  The  filter  was  adjusted  for  a  center  frequency  of  20  and  a 
bandwidth  of  3,  The  lower  frequency  did  not  cause  the  servo  to  chatter. 

The  general  procedure  used  in  the  simulation  for  Cases  t  ~  20,  101.9, 
and  200  was  as  follows: 

1.  The  computer  was  set  up  for  condition  3  of  the  particular 

flight  case. 

2.  Static  and  dynamic  checks  were  conducted. 

3.  The  Reference  level  and  KsA  were  set  at  the  values  used 
for  the  corresponding  aerodynamic  condition  and  dither  frequency  of  Flight 

Case  t  -  75. 

4.  The  values  of  Kg  as  read  from  the  computer  voltmeter  were 
recorded  for  aerodynamic  conditions  1,  2,  and  3  with  a  dither  frequency 
of  30  radians  per  second  and  for  condition  3  with  a  dither  frequency  of 

20  radians  per  second. 

5.  The  value  of  KsA  was  adjusted  to  give  a  suitable  rate  of 
change  of  Kg  for  the  particular  aerodynamic  condition  and  flight  case. 

The  Reference  level  was  then  adjusted  to  give  the  corresponding  value  of 
Kg  as  determined  in  Step  4  with  a  dither  frequency  of  30. 

6.  Recordings  were  made  at  the  dither  frequency  of  30  for 

aerodynamic  conditions  1  and  3.  During  the  run  for  each  condition,  Kg 
was  offset  by  a  step  voltage  put  into  the  integrating  servo  amplifier;  an 
actuator  signal,  ,  was  put  into  the  system  by  a  step  voltage  into  the 
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simulation  for  the  actuator;  a  wind  gust  was  simulated  by  a  step  voltage 
into  the  attack  angle,  o(  ,  amplifier;  and  a  dither  signal  failure  was 
simulated  by  turning  off  the  signal  generator. 

Step  6  was  also  conducted  for  Case  t  =  75. 

3. 3  Comparison  of  System  Selected  Gains. 

The  values  of  the  adaptive  gain.  Kg  9  which  were  recorded  by  the 
procedures  described  in  Section  3.2  during  the  analog  computer  simulation 
are  compared  in  this  Section  with  values  of  gain  determined  from  digital 
methods  to  establish  a  correspondnece  between  the  two  procedures.  A  Recomp 
II  Computer  was  used  for  the  digital  calculations.  The  damping  ratios  of 
the  short  period  modes  corresponding  to  the  recorded  adaptive  gains  are 
compared  to  establish  the  relative  adaptibility  of  the  system  under  vary¬ 
ing  aerodynamic  conditions  and  dither  frequencies  and  are  compared  with 
the  damping  ratios  corresponding  to  a  fixed  gain  system  to  establish  the 
benefit  of  the  adaptive  system  for  the  airframe  being  used  for  the  investi¬ 
gation,  The  analog  computer  simulation  was  conducted  for  aerodynamic 
conditions  1,  2,  and  3  as  described  in  Section  2.4.  The  digital  investiga¬ 
tion  included  aerodynamic  condition  4,  and  the  damping  ratios  associated 
with  the  gains  fcr  a  particular  Dither  Reference  level  for  this  aerodynamic 
condition  are  included  in  the  comparison  of  the  damping  ratios  of  the  ana¬ 
log  computer. 

The  determination  by  digital  methods  of  the  values  of  Kg  which  the 
system  would  select  was  accomplished  by  computing  the  amplitudes  of  the 
closed  loop  response  for  various  gams  and  for  frequencies  of  20  and  30 
radians  per  second.  Plots  of  the  amplitudes  of  the  closed  loop  response 
versus  for  a  dither  frequency  of  30  were  made  of  the  flight  cases  for 
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each  aerodynamic  condition.  For  aerodynamic  condition  3  plots  were  made 
for  dither  frequencies  of  30  and  20.  The  plots  are  shown  on  Figures  B-2 
(a)  through  B-2(f).  An  adaptive  gain  of  1  was  set  on  the  analog  computer 
for  each  aerodynamic  condition  of  Flight  Case  t  =  75  and  was  used  as  the 
standard  for  determining  the  values  of  Reference  level  and  adaptive  loop 
gain  pot  settings  which  were  maintained  constant  in  recording  the  adaptive 
gains  of  the  remaining  flight  cases.  The  adaptive  gain  value  of  1  for  Case 
t  =  75  was  also  used  as  the  standard  for  the  digital  method.  The  reference 
amplitude  of  the  system  was  determined  on  the  plots  by  the  intersection  of 
the  Kg  equal  to  1  line  with  the  Flight  Case  t  =  75  curve.  The  reference 
amplitude  eliminated  the  necessity  of  evaluating  the  components  of  the 
adaptive  loop  to  determine  a  Reference  Signal  level  and  was  used  to  deter¬ 
mine  the  values  of  gain  which  the  system  would  select  for  the  other  flight 
cases  by  the  intersection  of  the  reference  amplitude  and  the  particular 
flight  case  curve. 

The  comparison  of  the  adaptive  gains  determined  by  analog  and  digital 
methods  is  shown  in  Table  I  and  II  for  aerodynamic  conditions  land  2 
respectively  at  a  dither  frequency  of  30.  The  reference  amplitude  deter¬ 
mined  for  aerodynamic  condition  3  at  a  dither  frequency  of  30  was  too  large 
for  use  with  Flight  Cases  t  =  101.9  &  200  as  shown  on  Figures  B-2(c).  The 
dither  frequency  was  varied  in  the  digital  method  for  each  flight  case 
until  a  reference  amplitude  was  determined  which  resulted  in  the  gains  which 
were  selected  by  the  analog  computer.  The  dither  frequencies  required  to 
obtain  the  gains  are  shown  in  Table  III.  The  plots  of  closed  loop  response 
amplitude  versus  for  the  different  dither  frequencies  required  are  shown 
in  Figure  B-2(d).  The  third  bending  mode  was  excited  at  a  dither  frequency 
of  30  radians  per  second  and  the  filter  bandwidth  allowed  the  bending  signal 
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TABLE  I 


Gain  Comparison  for  Aerodynamic  Condition  1  at  Dither  Frequency  of  30 


j 

Flight  Case 

Gain  Selected  by 

Analog  Computer 

Gain  Determined  by 
Digital  Procedure 

t  =  20 

1.73 

1 

1.69 

t  =  75 

1.0 

1.0 

t  =  101.9 

0.54 

0.535 

t  =  200 

1 

1.15 

1.  17 

TABLE  II 


Gain  Comparison  for  Aerodynamic  Condition  2  at  Dither  Frequency  of  30 


Flight  Case 

Gain  Selected  by 

Analog  Computer 

Gain  Determined  by 
Digital  Procedure 

rr 

(I 

M 

O 

1.35 

1.3 

t  =  75 

1.0 

1.0 

t  =  101.9 

0.8 

0.72 

t  =  200 

3.4 

3.42 

3A 


to  pass  through  the  adaptive  loop  as  was  discussed  in  Section  3.2. 

The  dither  signal  source  could  not  be  accurately  set  to  within  0.628 
radians  per  second,  and  since  two  dither  frequencies  were  used  for  each 
flight  case,  the  signal  source  had  to  be  adjusted  for  each  case.  The 
variation  in  the  digital  frequencies  required  to  obtain  the  analog  gains 
is  approximately  within  the  accuracy  of  the  signal  source  settings.  The 
necessity  for  very  accurate  determination  and  setting  of  the  dither  fre¬ 
quency  and  Reference  Signal  level  when  bending  modes  are  excited  close  to 
the  dither  frequency  is  illustrated  by  the  difficulty  in  obtaining  a  compari¬ 
son  of  the  gains  for  the  digital  and  analog  methods*  When  bending  or  noise 
signals  pass  through  the  adaptive  loop,  there  is  a  possibility  of  the  system 
having  a  choice  of  selecting  either  of  two  adaptive  gains  which  are  relative¬ 
ly  close  in  magnitude.  The  selection  of  the  wrong  gain  might  cause  the 
system  to  be  unstable.  As  is  shown  in  Figure  B-2(d),  the  reference  ampli¬ 
tude  intersects  the  flight  cases  at  two  relatively  close  values  of  Kg 
The  comparison  of  the  gains  determined  for  aerodynamic  condition  3  at  a 
dither  frequency  of  20  is  shown  in  Table  IV.  The  difficulties  experienced 
with  a  dither  frequency  of  30  were  not  experienced  at  a  frequency  of  20 
where  bending  modes  are  not  excited.  The  plots  of  closed  loop  response 
versus  for  a  dither  frequency  of  20  are  shown  in  Figure  B-2(e), 

Root  loci  were  determined  for  the  various  flight  cases  and  aerodynamic 
conditions  using  a  Recomp  II  digital  computer.  The  damping  ratios  of  the 
short  period  modes  corresponding  to  the  gains  selected  by  the  analog  computer 
with  a  dither  frequency  of  30  and  to  the  gains  selected  by  digital  techniques 
for  aerodynamic  condition  4  with  a  dither  frequency  of  30  were  determined 
from  the  root  loci  and  are  compared  in  Table  V. 
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TABLE  III 


Frequencies  Required  by  Digital  Method  for  G"ins  Determined  by  Analog 
Method  for  Aerodynamic  Condition  3  at  Dither  Frequency  of  30 


| Flight  Case 
_ 

Gain  Selected  by 

Analog  Computer 

Frequency  Required 
by  Digital  Method 
in  Radians  per  Second 

o 

CM 

II 

U 

1.0 

28.5 

j  t  =  75 

1.0 

28.5 

* 

i  t  =  101.9 

j 

2.  1 

f 

30.0 

;  j 

t  =  200 

5.0 

N) 

00 

'-J 

TABLE  IV 

Gain  Comparison  for  Aerodynamic  Condition  3  at  Dither  Frequency  of  20 


Flight  Case 

Gain  Selected  by  Analog 
Computer 

Gain  Determined  by 
Digital  Procedure 

=  20 

1.35 

1.42 

t  -  75 

1.0 

1.0 

t  =  101.9 

j 

0.6 

j 

0.6 

t  =  200 

1 

1.5 

1.5 

1 
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TABLE  V 


Damping  Ratio  Comparison  for  Flights  Under  Various  Aerodynamic  Conditions 


1,1  . . . 

Flight  Case 

Damping  Ratios 

Damping  Ratios 

Damping  Ratios 

Damping  Ratios 

for 

for 

for 

for 

Aerodynamic 

Aerodynamic 

Ae rodynamic 

Aerodynamic 

Condition  1 

Condition  2 

Condition  3 

Condition  4 

t  =  20 

.82 

.78 

.725 

.685 

t  =  75 

.794 

.838 

.843 

.72 

t  =  101.9  i 

j 

.838 

1 

.87 

.912 

.875 

t  =  200 

! 

.89 

|  .  . .   .  .  . 

.924 

.93 

.89 

Total  1 

j  j 

i 

Var iat ion 

1 

in  Damping  \ 

1  i 

•  •  *■ 

\ 

.14 

.21 

.21 

R~t  io 

L 

Average 

| 

Damping 

i 

.792 

Ratio  for 

.836 

* 

.853 

.852 

Flight 

< 

t 

~  "  * .  — ~- 

t 

Variat ion 

i 

I  from  t  =  75 

i  + . 09 : 

i 

% 

1 

+  .086 

+  .087 

i  +.17 

J 

Reference 

1  -.o 

s 

4 

L 

-.058 

-.118 

1  -.035 

L - 

1 

4 -  — - 
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Flight  La t  75  as  jsed  as  the  reference  for  both  the  analog  and 
digital  procedures;  and  therefore,  the  maximum  variations  of  the  damping 
ratios  of  a  particular  flight  from  the  damping  ratio  of  cases  t  =  75  is 
probably  the  best  criterion  co  show  the  effect  of  bending  modes  on  the 
adaptibility  of  the  system.  The  results  of  the  comparison  of  the  fipst 
three  aerodynamic  conditions  shows  adaptibility  decreasing  as  bending  is 
included,  but  the  comparison  is  valid  only  for  a  dither  frequency  of  30 
and  must  be  checked  for  other  dither  frequencies.  The  investigation  as 
to  the  effect  on  the  results  of  a  comparison  when  the  dither  frequency  is 
varied  is  contained  in  Chapter  4.  A  different  method  was  used  in  aero¬ 
dynamic  condition  4  than  in  condition  3  to  remove  the  second  bending  mode. 

The  different  methods  resulted  in  a  large  variation  of  the  position  of 
the  zeros  associated  with  the  third  bending  mode  as  is  noted  by  examining 
the  transfer  functions  in  Appendix  A,  The  comparison  of  the  two  conditions 
in  Table  V  indicates  that  for  a  dither  frequency  of  30  the  frequency  track¬ 
ing  and  notch  filter  scheme  for  removing  the  second  bending  mode  results  in 
a  more  adaptive  system  than  the  bending  cancellation  scheme.  The  total  varia 
tion  of  the  damping  ratios  of  the  two  conditions  is  the  same  but  condition  4 
includes  a  slosh  mode,  has  the  larger  value  of  variation  from  the  reference 
value  of  the  t  =  75  case  towards  greater  damping,  and  has  an  average  damping 
ratio  during  the  flight  less  than  condition  3.  For  a  larger  average  damping 
ratio  during  a  flight,  a  smaller  variation  in  the  damping  ratios  would  be 
expected  because  a  larger  average  damping  ratio  means  larger  gains  are  used 
and  large  gains  force  the  damping  ratios  towards  the  constant  damping  ratio 
determined  by  the  position  of  the  Inverse  Model. 

The  damping  ratios  as  determined  from  the  root  loci  for  the  gains  select 
ed  by  the  analog  computer  for  aerodynamic  condition  3  with  dither  frequencies 
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of  30  ar.3  7f  .?  i-  iabl e  VI  with  the  damping  ratios  of  a  fixed 

gain  system  vith  the  c>  it  set  at  a  value  of  1.  The  gain  of  a  value  of  1 
was  selected  for  the  fixed  gain  system  in  order  to  establish  the  same 
reference  for  Case  t  -  75  that  was  used  for  the  Dither  System*  The  results 
of  the  comparison  show  that  the  fixed  gain  system  is  more  adaptive  than 
the  dither  system  at  a  frequency  of  30  radians  per  second  but  is  less 
adaptive  than  the  dither  system  at  a  frequency  of  20,  The  great  difference 
in  the  performance  of  the  dither  system  as  the  dither  frequency  was  varied 
led  to  the  investigation  of  obtaining  the  optimum  dither  frequencies  for 
the  different  aerodynamic  conditions.  This  investigation  is  contained  in 
Chapter  4. 

3.4  Reaction  of  the  System  to  Disturbances. 

Aerodynamic  conditions  1  and  3  were  used  in  the  investigation  of  the 
rate  at  which  the  adaptive  gain  was  adjusted  and  the  reaction  of  the  adap¬ 
tive  loop  to  wind  gusts,  actuator  movements,  and  loss  of  dither  signal. 

The  investigation  was  conducted  with  a  dither  frequency  of  30  radians  per 
second  so  that  the  evaulation  of  the  adaptive  loop  for  aerodynamic  condition 
3  would  be  performed  under  the  unfavorable  conditions  of  having  bending 
signals  included  with  the  filtered  dither  signal  in  the  adaptive  loop. 

Since  a  complete  analysis  of  the  adaptive  loop  was  not  performed,  the 
adaptive  loop  gain  was  adjusted  for  each  flight  case  with  the  following 
cons ide rat  ions : 

1.  To  produce  a  suitable  rate  of  change  of  the  adaptive  gain, 

2.  To  show  the  variation  in  the  loop  performance  for  different 
flight  cases  with  the  same  setting  of  the  servo  gain  pot, 

3.  To  use  several  values  of  for  the  flight  cases  with  one 

value  producing  an  overshoot  in  the  response. 
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TABLE  VI 


Damping  Ratio  Con-par  1  son  of  Dither  System  and  Fixed  Gain  System 


Flight  I 

Case 

Damping  Ratios 
for 

Dither  Frequency 
of  30 

Damping  Ratios 
for 

Dither  Frequency 
of  20 

Damping  Ratios 
for 

Fixed  Gain  System 

t  =  20 

.725 

.78 

.725 

t  =  75 

.843 

.838 

.843 

t  =  101.9 

.912 

.844 

.88 

t  =  200 

.93 

.9 

.882 

Total 
Variation 
in  Damping 
Rat  io 

.21 

.  12 

.  16 

Average 
Damping 
Ratio  for 
Flight 

.852 

.84 

.832 

Variat ion 
from 
t  =  75 
Reference 

1 

1  +.09 

j  -.118 

4.06 

-.06 

+  .039 
-.118 
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The  t  .  i  *  ■:  £  v  r  corresponding  flight  cases  under  the 

difftJL'c  i  it  tons  to  show  the  effect  of  bending  on  the 

loop  per  for  at  --:.* 

Record!:  ^  were  made  of  the  reaction  of  the  system  to  various  distur¬ 
bances  for  the  different  flight  cases  and  aerodynamic  conditions.  The  traces 
for  the  flight  cases  for  aerodynamic  condition  3  showing  the  actuator 
command  signal the  attack  angle,  C>(  ;  the  adaptive  gain  response,  Kg  ; 
the  signal  from  the  adaptive  loop  rectifier,  SAE ;  and  the  first  bending 
mode,  dp  are  contained  in  Figures  B-3(a)  through  B-3(d).  The  procedure 
used  in  obtaining  the  traces  is  contained  in  Steps  5  and  6  of  the  procedure 

used  to  simulate  Cases  t  -  20,  101*9,  and  200  in  Section  3*2. 

The  comparison  of  the  rates  of  change  of  Kg  when  Kg  is  offset  above 
and  below  the  undisturbed  magnitudes  selected  by  the  system  for  the  vari¬ 
ous  cases  and  conditions  is  shown  in  Table  VII*  Cases  t  -  20  and  75  have 
the  same  setting,  or  Kg  to  servo  input  ratio,  but  th,.  response  times 

vary  considerably*  The  variation  illustrates  the  fact  that  the  Kg  to 
servo  input  ratio  must  also  be  programmed  for  the  flight  if  constant  dynamic 
conditions  are  to  exist  in  the  system.  There  must  be  an  optimum  Reference 
Level,  Dither  Frequency,  and  Adaptive  Loop  Gain  for  optimum  operation  of  the 
dither  system.  The  program  required  to  maintain  the  optimum  value  of  the 
Kg  to  servo  input  ratio  can  be  placed  into  the  loop  in  the  form  of  a  shaped 
gain  curve  for  pot  as  is  explained  in  Section  3*5*  The  adaptive  loop 

gain  for  Case  t  =  101.9  resulted  in  an  overshoot  in  the  Kg  response  as  is 
shown  on  Figure  B-3(c)  The  response  time  for  this  case  is  the  fastest  as 

would  be  expected  since  the  Kg  responses  for  the  other  cases  are  ovv^  damped. 

As  bending  was  added  the  response  time  increased  except  for  Case  t  ~  101*9 
where  the  response  time  remained  constant.  In  order  to  form  some  basis 
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TABLE  VII 


Compar  son  u  katt j  uf  Change  of  the  Adaptive  Gain  for  Various  Loop  Gains 

and  Aerodynamic  Conditions 


Flight 

Case 

Undisturbed 

Values  of 

Time  (Seconds) 
per  Unity  Change 
in  K* 

Gain  Setting  in 
Adaptive  Loop 

ksa  k« 

Aerodynamic 

Condition 

1  3 

1  3 

j 

offset  aboo" 
t  =  ZU  j 

offset  belowj 
i 

!  -,o 
[ 
i 

3.4  22.0 

3.8  4.5 

.025  .025 

.025  .025 

.  i 

offset  above 

t  =  75 

offset  below 

1.0 

5.0  14.0 

5.0  15.0 

.025  .025 

.025  .025 

offset  above 
t  -  101.9 

offset  below 

2.1 

2.6  2.6 

4.0  4.0 

0.2  0.2 

0.2  0.2 

offset  above 

t  =  200 

offset  belo\ 

5.0 

6.6  11.0 

4.4  7.2 

0.0  0.1 

0.0  0.1 
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for  establishing  if  the  t  ii  .as  listed  m  Table  VII  are  sufficiently  fast  for 
the  system  to  function  properly,  the  time  intervals  between  flight  cases 
are  divided  by  the  gain  changes  required  between  flight  cases  to  estab¬ 
lish  figures  of  merit.  The  figures  of  merit  are  based  on  the  assumption 
that  is  changing  linearly  with  the  time  between  flight  cases  and  are 

listed  below: 

t  =  20  to  75  ££>seconds  per  unitychange  in  K^’ 

t  -  75  to  101.9  24.6  seconds  per  unity  change  in  Kg 

t  =  101.9  to  200  33.8  seconds  per  unity  change  in  Kg' 

Even  though  the  loop  was  not  optimized  the  times  listed  in  Table  VII  are 
well  within  the  values  established  as  figures  of  merit. 

The  comparison  of  the  reaction  of  the  adaptive  loop  to  actuator  signals 
for  the  various  flight  cases  and  aerodynamic  conditions  is  shown  in  Table 
VIII.  The  magnitude  of  K£  was  changed  in  all  cases  and  conditions  when  an 
actuator  signal  was  introduced  into  the  system,  but  the  adaptive  loop  re¬ 
turned  I  g  ,  while  the  signal  was  still  being  applied,  to  tho  original 
magnitudes  for  all  cases  and  conditions  in  relatively  small  time  intervals. 
The  ideal  system  would  result  in  no  change  in  when  an  actuator  signal 
was  experienced.  Excessive  per  cent  changes  in  occurred  for  Flight 

Case  t  =  20  for  aerodynamic  conditions  1  and  3  and  for  Flight  Case  t  ~ 

101.9  for  aerodynamic  condition  3.  If  the  adaptive  loop  had  an  optimized 
gain  with  a  properly  shaped  gain  pot  instead  of  the  linear  pot  used  in 

the  simulation,  these  large  changes  in  Kg  would  be  greatly  reduced.  One 
of  the  considerations  that  must  be  used  in  arriving  at  the  proper  shape  for 
the  gain  curve  of  the  pot  is  to  have  a  small  ratio  of  to  servo  input  for 
small  values  of  %  The  undisturbed  magnitudes  of  Kg  for  Case  t  a  20  was 
1  and  for  t  a  101.9  was  2.1  so  that  the  shaped  pot  would  aid  these  cases 
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table  vru 


Coir, par  i  sor.  u  f  the  Reaction  of  the  Adaptive  Loop  to 
Actuator  Signals 


Flight  Case 
Aerodynamic 
Condit ion 

3 

<Sc  Step 

Input  in 
Degrees 

Change  in 
Angle  of 
Attack,  0(y 
in  Degrees 

Change  In 

u 

Per  Cent 
Change  in 

Kj- 

• 

Time(Seconds) 
to  return  to 
Undisturbed 
Value  of 

Ga  in 

Setting 
in  Adaptive 
Loop  Ksa 

t  =  20 

-0.6 

1.0 

-0.2 

20.0 

4.6 

.025 

+0.6 

1.0 

-0.2 

20.0 

4.8 

.025 

t  =  75 

-0.4 

0.6 

-0.08 

8.0 

5.0 

.025 

1 

+0.4 

0.6 

+0. 16 

16.0 

20.1 

.025 

t  =  101.9  ' 

-0.15 

0.6 

+  0.8 

33.0 

4.8 

0.2 

+0. 15 

0.6 

-0.8 

33.0 

3.0 

0.2 

t  =  200 

-0.62 

0.5 

-0.4 

7.1 

9.0 

0.1 

+0.62 

0.5 

+0.4 

7.1 

8.5 

0.1 

Flight  Case 

Aerodynamic 

Condition  1 

t  =  20 

-0.85 

1.1 

-0.2 

20.0 

3.0 

.025 

+0.85 

1.1 

+0.2 

20.0 

3.0 

.025 

t  =  75  | 

-0.7 

0.7 

-0.11 

10.0 

3.6 

.025 

+0.7  ; 

0.7 

+0.11 

10.0 

3.8 

.025 

t  =  101.9 

-3.2 

1.0 

+0.16 

7.2 

3.0 

0.2 

+0.2 

1.0 

-0.08 

3.6 

2.0 

0.2 

t  =  200 

-0.5 

0.45 

-0.2 

3.8 

8.0 

0.1 

|  1 

+0.5 

! 

0.45 

+0.2 

3.8 

7.8 

0.1 

in  p.eventUig  lacgc  m  k g  when  disturbances  were  experienced. 

The  same  magnitude  actuator  signals  for  the  bending  aerodynamic  condition 
produced  approximately  the  same  change  in  the  angle  of  attack  as  for  the 
no  bending  aerodynamic  condition*  but  the  bending  condition  resulted  in 
greater  changes  in  and  longer  time  intervals  to  return  to  conditions 

that  existed  before  the  disturbance  than  the  no  bending  condition. 

The  comparison  of  the  reaction  of  the  adaptive  loop  to  wind  gusts  for 
the  various  flight  cases  and  aerodynamic  conditions  is  shown  in  Table  IX, 

A  one  degree  change  in  the  angle  of  attack  corresponds  to  wind  velocities 
greater  than  would  be  expected  by  the  missile  during  the  flight.  The  per 
cent  change  in  K^*  for  all  flight  cases  and  conditions  is  small,  within  11 
per  cent.  The  values  of  Kg  are  returned  to  the  undisturbed  magnitudes 
while  the  gust  is  still  being  experienced  for  all  flight  cases  and  condi¬ 
tions  in  small  intervals  of  time.  The  per  cent  change  in  and  the 

time  interval  to  return  to  conditions  that  existed  before  the  disturbance 
increase  as  the  airframe  becomes  more  flexible. 

The  Reference  signal  is  set  so  that  if  the  dither  signal  source  fails 
the  value  of  K£  is  increased  to  insure  that  the  system  will  not  become 
unstable.  The  comparison  of  the  rates  of  change  of  for  the  various 

flight  cases  and  aerodynamic  conditions  when  a  dither  failure  occurs  is 
shown  in  Table  X,  The  rate  of  increase  of  decreases  as  the  airframe 

becomes  more  flexible. 

The  results  of  the  investigation  contained  in  this  Chapter  show  that 
for  the  computer  simulation  of  the  dither  system  the  performance  of  the 
adaptive  loop  deteriorated  when  the  airframe  became  flexible.  The  adaptive 
gain  is  varied  by  actuator  signals  and  wind  gusts  while  the  ideal  system 
would  have  no  variation  of  gain  when  experiencing  disturbances.  The 
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.niiLE  IX 


^  iiparisovi  the  Reaction  of  the  Adaptive  Loop  to  Wind  Gusts 


|  Flight  Case 
Aerodynamic 
Condition  3 

Magnitude  of 
Change  in 

Angle  of 
Attack  in 
Degrees 

Change  in 

Per  Cent 
Change  in 

Time (seconds) 
to  Return  to 
Original 

Gain 

Setting 
in  Adan- 
tive  Loop 

t  *  20 

-4.0 

0.0 

0.0 

0 

.025 

+4.0 

C  .0 

0.0 

0 

.025 

t  =  75 

-1.0 

+0.04 

4.0 

6.4 

.025 

+  1.0 

-0.1 

10.0 

6.4 

.025 

t  «  101.9 

-1.0 

-0.4 

1.7 

2.4 

0.2 

+  1.0 

-0.4 

1.7 

2.4 

0.2 

t  =  200 

-3.0 

-0.6 

10.7 

11.0 

0.1 

+  3.0 

-0.6 

10.7 

14.6 

0.1 

-  Flight  Case 

Aerodynamic 

Condition  1 

t  =  20 

-1.0 

0.0 

0.0 

0 

.025 

+  1.0 

0.0 

0.0 

0 

.025 

t  =  75 

-1.0  ' 

+0.02 

1.8  ! 

4.0 

.025 

+  1.0 

-0.09 

8.2 

4.4 

.025 

t  =  101.9 

-1.0 

0.0 

0.0 

0 

0.2 

+  1.0 

+0. 1 

4.5 

1.2 

0.2 

t  =  200 

-3.0 

|  0.0 

0.0 

0 

0.  1 

+  3.0 

|  0.0 

0.0 

0 

0.1 

i  1 

i 

<  _ _ _  ___ 
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iABLE  X 


Comparison  of  Rates  of  Change  of  the  Adaptive  Gain  for  Dither  Failures 


Flight  Case 

Aerodynamic  Condition  3 

K£  Increase  Per  Second 

t  =  20 

1.2 

t  =  75 

1.12 

t  =  101.9 

5.6 

t  =  200 

1.8 

Flight  Case 

Aerodynamic  Condition  1 

t  =  20 

1.2 

t  =  75 

3.8 

t  =  101.9 

6.2 

t  =  200 

3.2 

variations  of  gain  for  wind  gusts  were  small  and  were  corrected  during 


the  gusts  in  a  relatively  short  period  of  time.  The  variations  of  gain 
for  actuator  signals  were  large  for  two  flight  cases  and  small  for  the 
remaining  flight  cases.  The  variations  of  gain  were  corrected  during 
the  signal  in  relatively  shor*  periods  of  time.  The  rates  of  change  of 
adaptive  gain  were  satisfactoi  ly  based  on  the  criterion  of  a  constant 
rate  of  change  of  gain  required  between  flight  cases.  By  optimizing 
the  gain  of  the  self  adaptive  loop  and  properly  shaping  the  adaptive 
gain  pot  the  rates  of  change  of  gain  could  be  increased*  the  variation 
in  the  adaptive  gain  due  to  actuator  and  wind  gusts  could  be  greatly 
reduced,  and  the  time  required  to  correct  for  disturbances  could  be 
reduced . 

3.5  Self  Adaptive  Loop, 

The  importance  of  optimizing  the  gain  of  the  adaptive  loop  was  shown 
in  Section  3.4.  For  small  perturbations  the  loop  is  linearized  in  order 
that  a  qualitative  investigation  may  be  conducted  (8).  The  linearized 
loop  with  the  Reference  value  as  the  input  is  shown  in  Figure  3.1,  The 
Reference  value  is  compared  with  the  rectified  filtered  dither  signal, 
and  an  error  signal  is  formed.  The  error  signal  drives  the  integrating 
servo  with  the  output  shaft  position  indicated  by  (j)  „  The  adaptive  gain 


which  is  the  slope  of  the  curve  of 


versus  servo  shaft  position.  The  gain,  ,  is  then  multiplied  by 


the  Closed  Loop  Response  Amplitude  versus  K ^  curve  at  the  Reference 
value  for  the  particular  flight  under  consideration.  The  resulting  signal, 


0,  is  then  fed  back  through  the  filter.  The  characteristic  equation  of  the 
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INTEGRATING  '  ADAPTIVE  VEHICLE 

SERVO  1  GAIN  POT  DYNAMICS 


EPPECT I VE  DYN AMI C  S 
OP  BANDPASS  PIETER 


i  Figure  3.1.  Linearized  Sell'  Adaptive  Loop. 


1.  The  root  locus  Is  shown  in  Figure  3.1. 


As  the  loop  gain  increases,  the  response  oscillates.  It  was  found  dur¬ 
ing  the  analog  computer  study  that  if  the  gain  of  the  adaptive  loop  was 
increased  beyond  a  certain  value  for  a  particular  flight  case  that  the 

response  would  oscillate  and  as  the  gain  was  reduced  the  Kg  response 
became  over  damped.  The  observations  of  the  study  agree  with  the  root 
locus . 

In  order  to  maintain  the  adaptive  loop  gain  constant  at  the  optimum 
value  that  is  determined,  the  product  of  the  Adaptive  Gain  Pot  transfer 
function  and  the  Vehicle  Dynamic  transfer  function  must  be  a  constant.  If 


at  the  value 


the  pot 


is  the  reciprocal  of 


of  K£  that  the  system  selects  for  the  particular  flight  case,  the  value 
of  Kv  will  be  a  constant  equal  to  1.  The  shape  of  the  gain  curve  for  the 
pot  for  aerodynamic  condition  3  with  a  dither  frequency  of  30  radians  per 
second  is  shown  in  Figure  3.2. 

In  an  actual  design  a  fixed  loop  gain  could  be  established  by  varying 
and  Kp  until  a  desirable  time  response  was  obtained  for  the  servo  out¬ 
put  shaft,  and  the  gain  of  the  pot  could  be  shaped  as  explained  above. 

If  the  performance  of  the  loop  was  not  satisfactory  as  to  the  rate  of  gain 
change,  the  change  in  gain  caused  by  disturbances,  or  the  time  required  to 
correct  the  gain  after  a  disturbance,  compromises  would  have  to  be  made. 

The  slope  of  the  gain  curve  could  be  varied  to  improve  the  performance. 

In  varying  the  slope,  however,  Kv  would  vary  and  the  adaptive  loop  gain  would 
no  longer  ba  a  constant  for  a  flight.  A  varying  Kv  would  vary  the  time 
response  of  the  servo  shaft,  and  the  performance  of  the  loop  would  be  af¬ 
fected.  A  different  dither  frequency  with  its  corresponding  Reference 
value  could  vary  the  shape  of  the  gain  curve  resulting  in  the  satisfactory 
performance  of  the  loop.  Optimum  setting  of  the  loop  gains  and  shape  of  the 
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for  the  k£  pot  could  be  determined  by  a  trial  and  error  procedure. 


gain  curvt 
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CHAPTER  4 


CLOSED  LOOP  RESPONSE  AMPLITUDE  AS  A  CRITERION  FOR  CONSTANT  DYNAMIC  RESPONSE 

4. 1  Introduct ion. 

The  object  of  this  Chapter  is  to  determine  for  the  dither  system  the 
accuracy  of  evaluating  the  dynamic  performance  of  a  variable  parameter  con¬ 
trol  system  -  a  flight  control  system  -  by  the  amplitude  of  the  closed  loop 
response.  A  dither  signal  is  used  as  an  input  to  the  system  and  is  varied 
by  the  dynamic  characteristics  of  the  system.  The  dynamic  characteristics 
are  then  evaluated  by  the  comparison  of  the  amplitude  of  the  dither  signal 
with  a  fixed  level  Reference  signal  which  represents  the  desired  dynamic 
characteristics.  The  optimum  dynamic  performance  of  the  system  under  in¬ 
vestigation  was  defined  in  Section  2,3  as  a  damping  ratio  of  the  short 
period  mode  of  0.7  with  a  minimum  variation  of  the  natural  frequency. 
Therefore,  the  accuracy  of  evaluating  the  dynamic  performance  by  the  ampli¬ 
tude  of  the  closed  loop  response  is  based  primarily  on  the  accuracy  of 
evaluating  the  damping  ratio  of  0.7. 

The  Reference  signal  values  for  dither  frequencies  of  20  and  30 
radians  per  second  were  not  determined  for  an  average  damping  ratio  of  0.7 
for  the  flight  for  the  computer  siim  letion  discussed  in  Chapter  3.  The 
Reference  values  used,  however,  resulted  in  approximately  equal  average 
damping  ratios  for  the  flights,  but  the  maximum  variation  of  the  damping 
ratios  was  twice  as  large  for  the  dither  frequency  of  30  as  for  the  frequency 
of  20.  The  maximum  variation  of  the  damping  ratios  of  a  flight  also  differed 
with  the  flexibility  of  the  airframe  and  with  the  consideration  of  variable 
fuel  load.  The  effect  of  varying  the  dither  frequency,  of  increasing  the 
flexibility  of  the  airframe,  and  of  varying  the  center  of  gravity  on  the 
accuracy  of  the  performance  evaluation  will  be  determined  in  this  Chapter. 

The  results  of  the  investigation  will  permit  the  selection  of  the  Reference 
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signal  vdlue  and  the  dither  frequency*  or  frequencies*  which  result  in 
the  optimum  variation  of  the  damping  ratio  and  natural  frequency  of  the 
short  period  mode  from  the  desired  system  dynamic  performance. 

The  investigation  was  conducted  principally  from  an  empirical  view¬ 
point  although  some  analytical  verification  is  included  in  Appendix  C. 
Aerodynamic  Conditions  1,  2,  3,  and  4  and  Flight  Cases  t  =  20,  75,  and 
101,9,  which  are  described  in  Section  2,4,  are  included  in  the  investigation. 
Flight  Case  t  =  200  was  not  included  because  the  sensors  are  located  in  a 
different  location  than  the  sensors  for  the  other  cases.  The  performance  of 
flights  evaluated  from  cases  that  are  measuring  the  reactions  of  the  systen 
at  different  positions  in  the  airframe  was  not  considered  a  true  indication 
of  the  capabilities  of  the  dither  system.  Therefore,  the  performance  of 
flights  are  evaluated  in  this  investigation  with  only  flight  cases  that 
are  measuring  reactions  at  the  same  airframe  position.  Application  of 
the  results  to  a  general  variable  parameter  system  is  discussed  in  Section 
4,4, 

4,2  Procedure. 

The  investigation  of  the  effect  of  varying  the  dither  frequency  on  the 
accuracy  of  the  closed  loop  response  amplitude  in  determining  the  damping 
ratio  of  0.7  was  conducted  using  Cises  t  =  20,  75,  and  101,9  as  representa¬ 
tive  of  a  typical  flight.  The  dither  frequency  was  varied  in  the  range 
from  0  to  40  radians  per  second.  The  upper  limit  of  40  radians  was  select¬ 
ed  because  available  actuators  will  not  respond  to  higher  frequencies.  For 
various  frequencies  in  the  dither  range  of  each  of  the  four  aerodynamic 
conditions  Reference  signal  values  were  determined  which  resulted  in  an 
average  damping  ratio  during  the  flight  of  0.7.  At  the  Reference  value  for 
a  particular  aerodynamic  condition  and  dither  frequency  the  maximum  varia¬ 
tion  of  the  damping  ratios  during  the  flight  was  determined.  Plots  of  the 
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maximum  damping  ratio  variation,/^  ,  versus  the  dither  frequency,  , 

were  then  made  for  the  four  aerodynamic  conditions.  The  procedure  employ¬ 
ed  in  determining  the  Reference  value  and  the  maximum  damping  ratio  varia¬ 
tion  will  be  described  by  the  use  of  a  typical  example.  A  dither  frequency 
of  22  radians  per  second  and  aerodynamic  condition  3  will  be  used  for  the 
example . 

Root  loci  for  the  three  flight  cases  were  determined  for  aerodynamic 
condition  3  using  a  Recomp  II  digital  computer.  The  damping  ratios,  £  , 

of  the  short  period  modes  associated  with  the  adaptive  gains,  Kg  ,  were 
determined  from  the  root  loci,  and  plots  of  S  versus  were  made  for  the 

flight  cases  as  shown  in  Figure  4.1.  The  amplitudes  of  the  condition  3 

closed  loop  transfer  functions  were  evaluated  at  the  dither  frequency  of 
22  radians  per  second  for  various  values  of  adaptive  gain.  The  plots  of 
closed  loop  response  amplitude  versus  adaptive  gain  for  the  flight  cases 
are  shown  in  Figure  4.2C  A  graph  of  the  damping  ratio  versus  flight  case 
was  constructed  as  shown  in  Figure  4.3.  The  procedure  used  in  establishing 

the  closed  loop  response  amplitude  which  would  result  in  an  average  damping 

ratio  during  the  flight  of  0.7  was  as  follows: 

1.  An  amplitude  for  the  closed  loop  response  was  selected,  and 

the  values  of  Kg  corresponding  to  the  amplitude  were  determined  for  the 
three  flight  cases  from  the  Amplitude  versus  plots. 

2.  The  values  of  the  damping  ratios  corresponding  to  the  values 

of  ,  which  resulted  from  the  selected  amplitude,  were  determined  from 

the  S  versus  plots. 

3.  The  values  of  damping  ratios  were  plotted  on  the  damping  ratio 
versus  flight  case  graph. 

4.  If  the  average  of  the  three  damping  ratios  corresponding  to 
the  same  closed  loop  response  amplitude  was  not  0.7,  a  different  response 

amplitude  was  selected.  Steps  1  through  4  were  repeated  until  the  average 
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X  Figure  4.2.  Closed  Loop  Response  Amplitude  vs  Adaptive  Gain  for  ^ 
Dither  Frequency  of  22  and  Aerodynamic  Condition  3  ill 
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damping  ratio  was  0.7 « 

When  the  closed  loop  response  amplitude  which  resulted  in  an  average 
damping  ratio  of  0.7  was  determined,  the  difference  between  the  maximum 
and  minimum  values  of  the  damping  ratios  for  the  three  flight  cases  was 
taken  as  the  maximum  variation  of  the  damping  ratio,  ,  during  the 

flight  for  a  dither  frequency  of  22  radians  per  second.  The  maximum 
variation  and  dither  frequency  were  plotted  on  the  Ac?  versus  VJ  graph, 
shown  in  Figure  4.4,  for  aerodynamic  condition  3.  The  procedure  was 
repeated  for  various  frequencies  for  the  various  aerodynamic  conditions 
to  construct  theA^  versus  2^/plots  for  the  dither  frequency  range. 

The  optimum  dynamic  performance  was  defined  as  a  damping  ratio  of  0.7 
with  a  minimum  variation  of  the  natural  frequency.  The  procedure  that  has 
been  described  determined  the  variation  of  the  damping  ratios  with  varying 
dither  frequencies  while  maintaining  an  average  ratio  of  0.7  during  a  flight. 
The  procedure  that  will  be  described  determined  the  variation  of  the  natural 
frequencies  of  the  short  period  modes  with  varying  dither  frequencies  while 
maintaining  the  same  average  damping  ratio  as  that  used  in  determining  the 
damping  ratio  variation.  The  reference  for  determining  the  damping  ratio 
variation  was  the  defined  optimum  performance  ratio  of  0.7.  The  reference 
for  determining  the  natural  frequency  variation  was  selected  as  the  average 
natural  frequency  of  the  short  period  modes  of  the  three  flight  cases  at  the 
damping  ratio  of  0.7.  The  procedure  employed  in  determining  the  natural 
frequency  variation  will  be  described  by  the  use  of  the  same  example  that 
was  used  for  the  damping  ratio  variation:  namely,  a  dither  frequency  of  22 
radians  per  second  and  aerc dynamic  condition  3. 

The  natural  frequencies  of  the  short  period  modes  at  a  damping  ratio 
of  0.7  were  determined  from  the  root  loci  for  the  three  flight  cases. 
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The  three  frequencies  were  averaged  to  determine  the  reference  natural 
frequency  for  condition  3.  The  associated  adaptive  gains*  K£  9  and 
natural  frequencies,  1U n,  of  the  short  period  modes  were  determined  from 
the  root  loci,  and  plots  of  ~U/c\  versus  were  made  for  the  flight  cases 

as  are  shown  in  Figure  4.5.  The  procedure  used  to  determine  the  frequency 
variation  was  as  follows: 

i 

1.  The  closed  loop  response  amplitude  which  resulted  in  an 
average  damping  ratio  of  0.7  for  the  flight  was  determined  from  the  damping 
ratio  versus  flight  case  graph  for  the  dither  frequency  of  22. 

2.  The  adaptive  gains  for  the  flight  cases  which  correspond  to 
the  amplitude  which  resulted  in  the  average  damping  ratio  of  0.7  were  deter¬ 
mined  from  the  Amplitude  versus  I<^  plots. 

3.  The  adaptive  gains  determined  from  the  Amplitude  versus 
plots  were  used  to  enter  the  1//^  versus  plots  to  determine  the  natural 
frequencies  of  the  three  flight  cases. 

4.  The  three  natural  frequencies  were  averaged.  The  reference 
natural  frequency  was  then  subtracted  from  the  averaged  frequency  to 
determine  the  natural  frequency  variation  for  the  dither  frequency  of  22 
radians  per  second. 

The  natural  frequency  variation,  ,  was  plotted  on  the 

versus  ~)jJ  graph  which  is  shown  in  Figure  4-6  for  aerodynamic  condition  3. 

The  procedure  was  repeated  for  various  frequencies  of  the  various  aerodynamic 
conditions  to  construct  the  AZ^versus  2t/  plot  s  for  the  dither  frequency  range. 

The  plots  of  versus  W  and  AliL.  versus  VJ  for  the  four  aerodynamic 
conditions  will  be  used  in  Section  4.3  to  determine  the  accuracy  of  the 
closed  loop  response  amplitude  as  a  criterion  for  constant  dynamic  response 
and  the  effect  on  the  accuracy  of  varying  the  dither  frequencies,  increasing 
the  flexibility  of  the  airframe,  and  varying  the  center  of  gravity. 


61 


cr* 


Figure  4.6.  Natural  1'requency  Variation  vs.  Ditner 
for  Aerodynamic  Conditions  1,  2,  3, 


frequency 
and  4. 


II 


M 


:rrt: 


63 


4 . 3  Results  of  the  Lmpiric-al  Investigation, 

One  of  the  basic  operations,  which  are  discuss  d  in  Section  2,1., 
of  an  adaptive  loop  is  the  continuous  measurement  of  the  system  dynamic 
performance.  The  measuring  method  must  have  a  negligible  effect  on  the 
system  response  signal.  In  the  dither  self  adaptive  flight  control  system 
the  measuring  method  involves  a  dither  signal*  and  the  important  system 
response  signal  is  the  short  period  mode  response*  The  possible  frequencies 
of  the  short  period  mode  response  are  never  greater  than  1*5  radians  per 
second  for  the  cases  and  aerodynamic  conditions  considered.  To  insure 
that  the  dither  signal  does  not  interfere  with  the  short  period  mode  re- 

ji 

sponse  the  lower  limit  for  the  dither  frequency  range  was  set  at  10  radians 
per  second.  The  upper  limit  was  previously  set  at  40  radians  per  second 
to  conform  to  the  characteristics  of  available  actuators. 

The  effect  of  varying  the  dither  frequency*  increasing  the  flexibility 
of  the  airframe,  and  varying  the  center  of  gravity  on  the  accuracy  of  evalu¬ 
ating  the  specified  damping  ratio  is  shown  in  Figure  4,4,  The  AS  versus 
7 U  plots  are  based  on  the  requirement  that  an  average  damping  ratio  of  0.7 
is  maintained  for  the  flight.  The  accuracy  with  which  the  damping  ratio  is 
evaluated  for  aerodynamic  conditions  1  and  2  is  independent  of  the  dither 
frequency  in  the  range  from  10  to  40  radians  per  second.  The  damping  ratio 
for  aerodynamic  condition  1,  which  is  the  rigid  airframe  condition,,  is  very 
accurately  evaluated  with  a  variation  of  0,06  between  the  high  and  low  values 
of  the  damping  ratios  during  the  flight.  As  the  airframe  becomes  flexible 
by  including  the  first  bending  mode,  aerodynamic  condition  2S  the  accuracy 
decreases.  The  variation  of  the  damping  ratios  for  condition  2  is  approxi¬ 
mately  0.15.  As  the  flexibility  of  the  airframe  increases  by  including  bend¬ 
ing  modes  1,  2,  and  3  but  eliminating  mode  2  by  a  bending  mode  cancellation 
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method,  aerodynamic  condition  3,  the  accuracy  of  evaluation  the  damping 
ratio  depends  on  the  dither  frequency,  A  difference  in  the  variation  of 
the  damping  ratios  for  dither  frequencies  of  20  and  30  radians  per  second 
was  determined  in  the  computer  simulation*  discussed  in  Chapter  3*  for 
this  aerodynamic  condition.  The  average  variation  of  the  damping  ratios 
for  the  dither  frequency  range  is  approximately  0*24  with  a  minimum  varia¬ 
tion  at  a  dither  frequency  of  23  radians  per  second  of  approximately  the 
same  magnitude  as  that  for  a  rigid  airframe,  namely*  0o06.  If  bending 
modes  1,  2,  and  3  with  bending  mode  2  eliminated  by  a  frequency  tracking 
and  notch  filter  method  and  a  varying  center  of  gravity  in  the  form  of  a 
slosh  mode  are  included  in  the  airframe  equations*  the  accuracy  of  evalua¬ 
ting  the  damping  ratio  is  independent  of  the  dither  frequency  from  12  to 
35  radians  per  second.  The  variation  of  the  damping  ratios  for  aerodynamic 
condition  4  in  the  12  to  35  frequency  range  is  approximately  0*435.  For 
dither  frequencies  greater  than  35  the  variation  decreases  and  is  0.375  at 
a  frequency  of  40  radians  per  second.  The  different  methods  used  in  eliminat¬ 
ing  the  second  bending  mode  in  aerodynamic  conditions  3  and  4  resulted  in 
large  differences  in  the  values  of  the  real  parts  of  the  zeros  associated 
with  the  third  bending  mode  of  the  airframe  transfer  function  as  is  discussed 
in  Section  3.2.  The  zero  for  the  third  bending  mode  for  condition  3  is  close 
to  the  imaginary  axis  while  the  zero  for  condition  4  is  not. 

In  summary,  the  frequency  in  the  dither  range  has  no  appreciable  affect 
on  the  accuracy  of  evaluating  the  damping  ratio  for  aerodynamic  conditions 
1,  2,  and  4.  The  accuracy  depends  greatly  on  the  choice  of  the  dither  fre¬ 
quency  for  condition  3  with  the  magnitude  of  the  variation  varying  from  the 
small  values  obtained  in  condition  1  to  the  large  values  obtained  in  condi¬ 
tion  4.  The  accuracy  decreases  as  the  airframe  becomes  flexible  and  as  a 


65 


varying  center  of  gravity  is  Included.  The  large  variations  of  condi¬ 
tion  4  could  be  reduced  by  using  an  average  damping  ratio  of  greater 
than  0.7  for  the  flight  because  for  large  damping  ratios  the  amount  of 
change  of  gain  per  unit  change  in  damping  ratio  is  greater  than  for  small 
damping  ratios  as  is  shown  in  Figure  4.1. 

The  effect  of  varying  the  dither  frequency^  increasing  the  flexibility 
of  the  airframe,  and  varying  the  center  of  gravity  on  the  natural  fre¬ 
quency  variation  is  shown  in  Figure  4.6.  The  variation  of  the  natural 
frequencies  for  aerodynamic  conditions  1  and  2  is  independent  of  the  dither 
frequency  in  the  range  from  10  to  40  radians  per  second.  The  variation  for 
condition  1  is  approximately  0.034  radians  per  second  and  for  condition  2 
is  approximately  0.123.  The  reference  natural  frequency  for  condition  1 
was  1.355  radians  per  second  and  for  condition  2  was  1.185.  The  variation 
of  the  natural  frequency  depends  on  the  dither  frequency  for  aerodynamic 
conditions  3  and  4,  The  average  variation  of  the  natural  frequency  for 
condition  3  for  dither  frequencies  in  the  dither  range  is  approximately 
0.068  radians  per  second.  At  dither  frequencies  of  approximately  23  and 
24  radians  per  second  there  are  no  variations.  The  average  variation  of 
the  natural  frequency  for  condition  4  is  approximately  0.06  radians  per 
second  with  a  minimum  variation  of  0.028  at  a  dither  frequency  of  40  radians 
per  second.  The  reference  natural  frequency  for  condition  3  was  1.225 
radians  per  second  and  for  condition  4  was  1.51. 

In  summary,  the  variation  of  the  natural  frequency  does  not  depend  on 
the  dither  frequency  for  aerodynamic  conditions  land  2  but  does  depend  on 
the  dither  frequency  for  conditions  3  and  4.  The  variations  were  small  for 
all  conditions  with  10  per  cent  of  the  corresponding  reference  value  being 
the  maximum  variation  determined.  The  effect  of  the  flexibility  of  the 
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airframe  and  of  a  varying  center  of  gravity  on  the  variations  of  the 

natural  frequency  is  not  clearly  indicated.  The  magnitudes  of  the  average 
variations  of  the  natural  frequencies  for  the  aerodynamic  conditions  varied 
from  large  to  small  in  the  following  order:  condition  3,  4,  and  1. 

The  investigation  of  the  variations  of  the  damping  ratios  and  natural 
frequencies  was  conducted  with  the  amplitude  of  the.  closed  loop  response 
selected  to  result  in  an  average  damping  ratio  daring  the  flight  of  0*7. 

Since  the  defined  optimum  performance  of  the  system  was  a  damping  ratio  of 
0.7  with  a  minimum  variation  of  the  natural  frequency ,  the  optimum  dither 
frequency,  or  frequencies,  for  each  of  the  aerodynamic  conditions  was 
determined  on  the  basis  of  minimizing  the  variation  in  damping  ratio  and 
natural  frequency.  The  error  criterion  used  to  determine  the  minimization 
of  the  variation  in  damping  ratio  and  natural  frequency  was 
Plots  of  AS  versus  are  shown  in  Figures  4.7  through  4.10  for  four 

aerodynamic  conditions.  The  optimum  dither  frequencies  selected  for  the 
aerodynamic  conditions  are  shown  in  Table  XI  with  the  corresponding  error 
criterion  value,  the  variation  in  the  damping  ratio*  the  variation  in  the 
natural  frequency,  and  the  Reference  signal  value.  Two  dither  frequencies 
with  different  error  criterion  values  are  listed  for  condition  4,  The  dither 
frequency  of  40  has  the  lower  criterion  value,  but  the  Amplitude  versus 
curve  for  that  frequency  has  a  very  small  slope  at  the  Reference  value  which 
would  require  a  very  accurate  setting  of  the  Reference  signal.  The  dither 
frequency  of  36  has  a  larger  slope  and  would  be  the  preferred  frequency  if 
the  system  was  actually  mechanized. 

The  Reference  signal  values  listed  in  Table  XI  are  equal  to  the  closed 
loop  response  amplitudes  required  for  an  average  damping  ratio  of  0.7  for 
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Optimum  Cither  Frequencies  and  R.  ier-rc-  Values 


Aerodynamic 
Condit ion 

Opt imum 

Dither 

Frequency 

or 

Frequencies 
(Rad/Sec0 ) 

Error 

Criterion 

Value 

Variat ion 
iv 

Ratio 

Variat ion 
in 

Natural 
Frequency 
{Rad/Sec. ) 

Reference 

Signal 

Value 

Condition  1 

10 

0.56 

26 

0,068 

0.06 

0.032 

0.485 

40 

0.42 

Condition  2 

10 

0.835 

12 

0,19 

C.  15 

0,118 

0.79 

26 

0.7 

Condition  3 

22 

0, 105 

0,095 

0.047 

0.43 

Condition  4 

40 

0.375 

0,375 

0.027 

0.993 

36 

0.405 

0,405 

0,028 

0.975 
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the  fli.hf  The  ^ « t  '\t  1  s  c  rreh.t  if  th*  .  jt  t*  r  /<  ..  alue  of 

the  dither  input  signal  nd  the  gains  of  tb  t.i  ■  v.  -nd  rectifier  in  the 
adaptive  loop  are  equal  to  unity.  If  the  Valdes  art  not  unity,  the 
Reference  value  listed  must  be  corrected  by  being  multiplied  by  the  pro¬ 
duct  of  the  half  cycle  average  of  the  filter,  and  the  gain  of  the  rectifier. 

4 . 4  Application  of  the  Results  to  a  General  Case,, 

The  object  of  this  section  is  to  formulate  a  procedure  from  the  empiri¬ 
cal  results  discussed  in  Section  4.3  and  from  the  analytical  verification 
contained  in  Appendix  C  to  facilitate  the  rapid  selection  of  the  optimum 
dither  frequency.  The  selection  will  be  based  >n  minimizing  the  damping 
ratio  variation.  Once  the  dither  frequency  or  area  of  frequencies  is 
selected,  the  procedures  of  Section  402  can  be  used  for  the  frequency  to 
determine  the  Reference  value  and  the  magnitude  of  the  damping  ratio  varia¬ 
tion  for  the  average  damping  ratio  desired  for  the  flight. 

The  dither  range  is  established  by  selecting  the  lover  frequency  limit 
five  or  more  multiples  greater  than  the  maximum  short  period  frequency  to 
eliminate  coupling.  The  maximum  frequency  of  the  short  period  mode  can  be 
approximated  by  the  square  root  of  the  square  of  the  real  and  imaginary 
parts  of  the  zeros  of  the  Inverse  Model.  The  upper  limit  of  the  range  is 
fixed  by  the  actuator  requirements  at  40  radians  pAr  second. 

The  assumption  is  made  that  the  control  system  will  be  investigated  by 
a  point  study  of  the  flight  trajectory.  From,  the  root  loci  shown  in  Figure 
4.11,  it  is  noted  that  the  open  loop  poles  associated  with  the  bending  modes 
remain  relatively  stationary  during  a  flight  but  that  the  pen  loop  zeros 
associated  with  bending  modes  move  during  the  flight  from  smaller  to  larger 
values.  If  slosh  modes  are  temporarily  Ignored,  chr  effect  that  the  varying 
values  of  the  bending  zeros  have  on  the  shore  period  mode  is  to  vary  the 
shape  of  the  short  period  characteristic.  The  bending  zer-s  must  be  in  the 
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Figure  4.11.  Root  Loci  Lor  Aerodynamic  Condition  3. 
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on  the  shc.-p-  Ihe  flight  jsc*  with  th  h  -.e  wbth  the  smallest 

imaginary  part**  which  is  the  earliest  flight  ...  ?nsid**red  after  launch, 
results  in  a  short  period  character  1st  1 :  whx.h  h  s  a  larger  adapative  gain  at 
a  particular  damping  ratio  than  any  of  the  other  cases*  The  open  loop  poles 
and  zeros  associated  with  the  short  period  appx  -•>  *  m;  r  i  on  nls  r-  vary  during 
the  flight,  but  it  is  assumed  that  the  variations  an1  negligible  compared 
to  the  variations  of  the  bending  zeros.  The  values  of  adaptive  gain  at  a 
particular  damping  ratio  decreases  as  the  imaginary  part  of  the  bending 
zero  increases.  The  adaptive  gains  for  a  particular  damping  ratio  are, 
therefore,  inversely  proportional  to  the  time  in  flight 0  The  dither  system 
maintains  the  amplitude  of  the  closed  loop  response  constant  daring  a  flight. 
A  constant  amplitude  results  when  the  product  of  the  adaptive  gain  and  the 
dynamic  gain  remains  constant.  The  dynamic  gain  is  the  gain  associated  with 
the  vehicle  dynamic  characteristics,,  In  order  to  have  a  minimum  variation  in 
the  damping  ratios,  the  product  of  the  dynamic  gain  and  the  adaptive  gain 
corresponding  to  the  desired  damping  ratio  must  be  approximately  constant 
during  the  flight.  Since  the  adaptive  gain  varies  inversely  with  the  time 
in  flight,  the  dynamic  gain  must  vary  directly  with  the  time  in  flight  for 
the  product  to  be  a  constant*  If  the  dither  frequency  is  selected  close 
to  a  bending  zero  of  the  initial  flight  ca sev  the  dynamic  gains  will  vary 
directly  with  the  time  in  flight.  The  herding  zero  must  have  a  value  with 
a  small  real  part  and  an  imaginary  part  in  the  dither  range. 

If  the  value  of  the  bending  zero  has  an  imaginary  part  in  the  dither 
range  but  a  large  real  part*  the  selection  of  the  dither  frequency  for 
optimum  damping  ratio  variation  is  independent  of  the  position  of  the  zero. 
The  value  of  the  dither  frequency  does  not  appreciably  affect  the  dynamic 

;s 


i*  1  .  •  »xir.a?t!v  •  sr.  •  r  d  mpi  g  xp  ;ed*  During 

the  fligh-  the  bending  ™ode  s  do  result  in  \o  Iti  ns  nr  :  ne  adaptive  gains 
for  the  same  damping  ratio,  Since  the  vai.,  t  r  dinner  frequency  can 
not  result  in  dynamic  gains  varying  inversely  t  the  short  period  gains* 
large  variations  in  the  damping  ratio  may  result  *  The  requirement  of  deter¬ 
mining  when  the  real  part  of  the  value  of  the  bending  zero  is  large  or  small 
can  be  eliminated  by  always  selecting  the  valup.  of  the  dither  frequency 
close  to  the  imaginary  part  of  the  bending  zero  of  the  initial  flight  case. 

If  bending  zeros  do  not  exist  In  the  dither  range ^  the  optimum  damp¬ 
ing  ratio  variation  is  independent  of  the  selection  of  the  dither  frequency. 
The  value  of  the  dither  frequency  can  not  vary  the  dynamic  gain  to  result 
in  adaptive  gains  which  correspond  to  the  values  required  for  constant  short 
period  damping  ratios.  The  variations  in  the  adaptive  gains  of  the  short 
period  modes  for  the  same  damping  ratio  increase  as  the  number  of  bending 
modes  increase;  and  therefore*  the  magnitude  of  the  damping  ratio  variation 
would  be  expected  to  increase  as  the  number  of  bending  modes  increase. 

Slosh  modes  may  have  an  important  effect  on  the  dynamic  performance 
of  the  system  for  some  airframes*  Slosh  modes  would  have  a  negligible 
effect  during  the  initial  and  final  phases  cf  the  flight  because  the  fuel 
tanks  are  rhen  approximately  full  or  empty.  The  slosh  modes  would  affect 
the  middle  phase  of  the  flight.  If  the  slosh  mode  frequencies  are  close 
to  the  short  period  frequencies y  the  shape  of  the  short  period  characteristic 
could  be  appreciably  varied  for  the  middle  flight  cases0  The  shape  of  the 
characteristic  would  not  be  affected  for  the  initial  and  final  flight  cases. 
The  largest  adaptive  gain  for  a  particular  damping  ratio  would  be  for  the 
middle  flight  case  which  is  most  affected  by  the  slosh  modes.  The  dynamic 
gain  would  have  to  be  reduced  by  the  selection  of  a.  dither  frequency  for 
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of  the  imaginary  parts  of  the  bending  aero  of  th.  *  initial  case  and  the  bend- 

ing  zero  of  the  middle  case  to  result  i*  a  small  damping  ratio  variation 

throughout  the  entire  flight e  The  exact  frequency  in  the  area  would  have  to 
be  determined  by  the  procedure  in  Section  402* 

In  summary,  the  selection  of  the  optimum  dither  frequency  is  as  follows! 

1.  If  no  bending  zeros  exist  in  the  dither  range,,  the  value  of  the 
dither  frequency  is  not  restricted* 

2.  If  there  are  no  important  slosh  modes  in  the  system  but  if 

there  are  bending  zeros  in  the  dither  range »  the  optimum  frequency  is  close 
to  the  imaginary  part  of  a  zero  associated  with  the  Initial  flight  case* 

3.  If  there  are  important  slosh  modes  in  the  system  and  there  are 

bending  zeros  in  the  dither  range s  the  optimum  frequency  area  is  between  the 
value  of  the  imaginary  part  of  an  initial  flight  case  zero  and  the  imaginary 
part  of  the  zero  of  the  middle  case  which  is  most  affected  by  the  slosh  modes* 

The  dither  frequency  must  not  assume  a  value  that  is  close  to  the  value 
of  the  imaginary  part  of  a  pole  which  is  adjacent  to  the  imaginary  axis  be- 
cause  the  adaptive  gain  change  per  unit  amplitude  of  the  closed  loop  re-* 
sponse  is  large,  A  large  gain  to  amplitude  ratio  requires  a  very  accurate 
setting  of  the  Reference  value* 
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CHAPTER 


CONCLUSIONS 

The  conclusions  from  the  results  of  the  investigation  are  as  follows: 

1.  The  dither  adaptive  loop  employs  a  single  parameter  to  control 
basically  a  single  dynamic  characteristic  of  the  system,,  The  parameter  of 
the  loop  is  the  closed  loop  response  amplitude  and  the  characteristic  con¬ 
trolled  is  the  damping  ratio, 

2.  The  components  of  the  loop  are  simple  to  mechanize^  but  a  prior 
knowledge  of  the  airframe  dynamics  are  required  to  select  the  dither  fre¬ 
quency,  Reference  value,  and  the  shape  of  the  gain  curve  for  the  adaptive 
gain  pot, 

3.  The  adaptive  loop  will  adjust  the  adaptive  gain  to  partially 
correct  for  the  varying  dynamic  characteristics  of  a  flexible  airframe  if 
the  dither  frequency  has  no  coupling  with  a  bending  mode,  The  accuracy 

of  the  correction  in  maintaining  a  constant  damping  ratio  is  the  inherent 
accuracy  of  the  Amplitude  Closed  Loop  Response  Criterion, 

4.  The  adaptive  loop  makes  satisfactory  corrections  of  the  adap¬ 
tive  gain  for  disturbances.  A  further  investigation  into  the  gain  curve 
shape  for  the  adaptive  gain  pot  and  into  the  optimum  loop  gain  is  required 
to  verify  that  an  actuator  signal  will  not  cause  an  unsatisfactory  variation 
in  the  adaptive  gain, 

5.  The  Amplitude  Closed  Loop  Response  Criterion  does  not  result 
in  a  constant  damping  ratio.  The  variation  of  the  damping  ratios  increases 
as  the  desired  average  damping  ratio  decreases*,  as  the  airframe  flexibility 
increases,  and  as  the  variation  of  the  center  of  gravity  increases.  If 
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frequencies  in  the  dither  range  greatly  affect  the  amplitude  of  the  vehicle 
dynamic  response,  the  variations  of  the  damping  ratios  can  be  appreciably 
reduced  by  the  proper  selection  of  the  dither  frequency*, 

6*  When  the  damping  ratio  variation  varies  with  the  dither  fre¬ 
quency,  the  minimum  natural  frequency  variation  occurs  at  approximately  the 
same  dither  frequency  as  the  minimum  damping  ratio  variation,, 
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APPENDIX  A 


The  leiodyiuirnir  equations  for  the  missile  were  obtained  from  the 
Autonetics  Division  of  North  American  Aviation,  The  equations  are  two 
degree  of  freedom  sh  rt  period  mode  approximations  with  three  bending 
modes  and  one  fuel  slosh,  mode  included. 
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I&dl  b  7)6^3  bfeds  -hL(=d.3-bM6ATl-h/(>AT-b06>Ai.6XJrPi>Ai.bbQ6F)T 


iuotfr/W  +67^4- C-}@1  O/di  rf?a//  -bFjJ)  T  G 7  dA-b  H?  bT? cA  4“ 

J7J5  b  Ryclz  f  i.7  q/s  4-  M7  Al  4*  /A  AioX  4-  Q?Al qX  b  FyAlH  b  (p7<^r 
A^/  =  Ai^bBsQbCsObD?d/-bE$<Ji’bFi?di-bG’sdA  4-  H$d&  -fTfriafoSJls* 
K?cta  bLsdib  MpAl  b  tV$Aic>xb  0%  A\lH  b  P^AlM  b  Q#At 


Gs-  =>  Q  b  /l /di  f  /L^Ja  b  /^3<d3 
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■  hcv 


cx  — 

©  - 

Jl  - 

rh  - 

Jj  - 

a .  >  1  e 

pitch 

first 

second 

third 

•u  ott. ■«  In  degrees 

angle  in  degrees 

bending  mode  in  feet 

.  bending  mode  in  feet 

bending  mode  in  feet 

©&  - 

pitch 

angle  as  measured  by 

position  gyro 

Ai  - 

first 

ignition  fuel  slosh  mode  in  feet 

AloX  - 

first 

liquid  oxygen  slosh  mode  in  feet 

At*  - 

first 

liquid  hydrogen  slosh 

mode  in  feet 

Jr  - 

thrust 

control  deflection  in  degrees. 

tj_20 

t-75 

t -101. 9 

t-200 

“ 

A, 

-.064 

-.0489 

-.0336 

-.0236 

B, 

1 

1 

1 

1 

C, 

-.158 

-.0226 

-.0103 

.0074 

D| 

0 

.00123 

.000229 

.00015 

E , 

0 

.0337 

.031 

-.061 

E, 

0 

-.000886 

-.000167 

-.00019 

G, 

0 

-.0461 

-.0154 

-.185 

Hi 

0 

.001108 

.000214 

.00014 

I| 

0 

.0940 

.0346 

-.047 

J, 

.0567 

.0289 

0 

0 

K| 

.058 

.0298 

.0422 

-.039 

L, 

.00497 

.00255 

.00362 

-.0027 

M, 

.222 

.0565 

.0378 

.02 

Afl 

0 

.0217 

.0159 

1.65 

3 .3 

0 

-.00156 

-.000519 

-.038 

0 

-.00204 

po 

-.000813 

-.055 

«,_2f) 

-  -  75 

tj_101 . 9 

t.^200 

0 

0515 

-.0386 

-31.92 

Ea 

0 

00168 

.000664 

,042 

F,a 

0 

.071 

.0552 

-12.73 

Ga 

0 

,00213 

-.000792 

-.024 

H* 

0 

-.163 

-.127 

-28.33 

Ift 

.00316 

.248 

0 

0 

-.0048 

-.00716 

.00193 

-3,8 

o*- 

-.000919 

-.00205 

-.00482 

.012 

.0485 

.0814 

.  157 

2,28 

a3 

0 

.401 

.213 

.32 

S3 

0 

-.0294 

-.00723 

-.012 

Cj 

-.163 

-.205 

-.175 

-.59 

D3 

-66.  5 

-67,3 

-66.9 

-767.46 

E3 

0 

.0348 

.00981 

.023 

F3 

0 

1.39 

.821 

16.34 

% 

0 

-.0482 

-.0125 

-.0013 

H3 

0 

-.3.24 

-1.83 

6.0 

T3 

.137 

-.0414 

0 

0 

J3 

.  131 

,298 

.634 

3.08 

K3 

-.00«01 

-.0182 

-.0389 

- .  128 

L3 

t  Pf  '• 

-  1  .  .  . 

-  i  _  ‘  ;v 

-1.79 

-1.86 

a* 

0 

-.288 

-.155 

-.39 

Bh 

0 

.0243 

.00591 

.009 

cH 

0 

.0348 

.00981 

.023 

D* 

0 

.708 

.384 

13.82 

E,,. 

-.43 

-.463 

-.439 

-.96 

F« 

-462 

-463 

-463 

-2169.63 
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r_7  5 

t-101.9 

t-200 

g4 

0 

.0459 

.0124 

.015 

V 

0 

2.85 

1.53 

12.66 

lH 

,21 

.373 

0 

0 

-.224 

-.51 

-1.086 

3.78 

K* 

-.0112 

-.0255 

.0544 

- .  146 

LV 

-1.57 

-1.78 

-1.79 

-2,4 

0 

.36 

.  199 

.28 

Br 

0 

-.0308 

-.00704 

-.0052 

cr 

0 

-.0482 

-.0125 

-.0013 

0 

- .  18 

-.43 

10.24 

Er 

0 

.0469 

.0.24 

.015 

0 

1.27 

.584 

12.7 

Gr 

-.832 

-.906 

-.852 

-1.49 

V 

-1730 

-1734 

-1732 

-5351.8 

V 

-.0182 

.476 

.18 

0 

Jr 

.0374 

.085 

.181 

3.52 

Kr 

.0288 

.0655 

-.14 

.018 

Lr 

-1.57 

-1.78 

-1.79 

-2.18 

X, 

-2.3 

-2,3 

-2.3 

0 

Xa 

3.4 

3.4 

3.4 

3.4 

Xj 

-6.95 

-6.95 

-6.95 

0 

The  magnitudes  of  the  remaining  coefficients  were  not  available  but  are 
included  in  the  transfer  functions. 
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Informat  ion 


1  ■  *x  t  i'-Uions  0  and  C\  will  be  shown. 

‘  0  %*  *  *  *  *  M 

L<  ,  't  ’'io  'if  ‘dti.»n  of  the  d,  s  ,  dj  ,  Ax  »  Au>* 
equation*,  is  contained  ir.  Reference  10* 


and  Aj.h 


From  Newton* s  second  law  of  motion. 


^  Fy  =  Ci'y  jrgr 

>*  T 

£L  -  3^-  i-MS7\  J, ) 

£M  =  +/(S7Xh) 

iff  -  Ay  +-£(ioxJ>) 


£F^  r  <£✓ 


(2) 


where  Fy  are  the  forces  in  the  x,  y,  z  directions  of  a 

right  hand  system  of  Cartesian  axes  fixed  in  space;  where  ^L*  ^  M, 
are  the  moments  about  the  x,  y,  and  z  axes;  where  h^,  ,  ,  and  h^  are 

the  moments  of  momentum  about  x^  ,  y(  ,  and  z^  axes  which  are  fixed  to 
the  airframe,  where  ^  (w  X  h  X  h  )  $  an  ^^(U/Xh^)  arise  from 

the  angular  velocity  'll/  of  the  Xj  ,  y^  ,  and  zj  axes  with  respect  to 
the  x,  y,  z  axes;  where  bl^c  ,  a^  ,  a^  ,  are  the  accelerations  in  the 
x,  y,  and  z  directions,  and  where  m  is  the  mass  of  the  airframe. 

The  absolute  acceleration  of  a  body  with  respect  to  x,  y*  z  axes  is 


aa/ra*  -  "F  ii/X  / 

where  V  is  the  instantaneous  linear  velocity  and  is  the  angular  velocity, 

U  +  V+-W 

M/Y  =  (Qv/-Rv)M  r  (mu- pw)m  t  (PV-ou)  JL 

where  P,  Q,  and  R  are  the  angular  velocities  along  the  xs  y,  and  z  axes 
respectively  and  U,  V,  and  W  are  the  linear  velocities  along  the  x,  y s  and 
z  axes.  The  components  of  the  acceleration  are  then 

~  U  tfyW-  R 1/ 
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/  ..  --  V  +  flU-PW 
rPV-QU 

Substituting  the  components  of  acceleration  into  Equations  1  and  consider¬ 
ing  the  mass  as  a  constant  the  equations  become 

CutQW-RV) 

(3)  £  Fy  =  **u(V  T  RU  -  pw ) 

F 2  -  a^-Cw  P  v "  cpt/  ,) 


The  components  of  the  moment  of  momentum  are  calculated  from  summing  the 
moments  of  the  velocity  vectors  about  each  axis  and  multiplying  by  the  mass 

-  PI*#.-  QX^^y  ~  RJ.yLsy'' 
hy  z  (j)J[yy  -  P?X y  ys  ~ PT 

where  I,?.^  r  ^(y'-f  y  A)  and  T-^y  ~  J  y  .y 

XY  and  XZ  are  assumed  planes  of  symmetry,  and  therefore 


(4) 


IsjL-y 


{n 


(5) 


(6) 


-  0.  Equations  A  become: 

h --PI 

h y  —  Cj)T  y  y 

Fy'  Flyy 

~  P 

M- 


-  £p  T/y/y 

wxl  =  0^  9  -  ^  (hf-R-h^fy  +  {hyP-^Q)J- 

Equations  2  with  Equations  5  and  6  substituted  become 


Since 
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(7) 


*  L  -  QR{l*.  ^ 

£m~  QJiyy  t  PRfof  ~~ 

<£  /\/  -  R  t  ^  ty(Tjj  —  X^6^'/)  • 

The  forces  indicated  by  the  left  terms  of  Equations  3  represent  the 

summation  of  aerodynamic,  thrust,  and  gravity  forces.  The  gravity  forces 

are  shown  below  and  their  derivation  can  be  found  in  Reference  9. 


(8) 


Coe-/  ^  0**  Ce<u  Qa 

O'* 


,) 


0  s&trts'  ^  ~  (s/TL'ps'  C&ds  Qo  Ccr4s  (fio )  sA&K*'  0 

*  ( o  )  (c*^  Q  —  £cr&-  ^/)f' 

£)n (c**'  ^ Ccr<is&tsterLS  (f)  f 

C U*ls  Qc  Ur^C^o  )(  Or**  &  (p) 

-  {^syrs'^-sO'U*/  &b )  i  Or^  W 40>\s  &  U>Oy  (p  +xul^  p)  ~b 

(  /in  -^'  tir^s  C^o  .j&Ohs'  s  &  Co^>  ~  Gru-  P d-tsrvCjPj  -f~ 


( 


/>  >1y-' 


■?' 


6^c>  X::))  (cey^  &  Cercu  (j) )  . 


where  m  is  the  mass  cf  the  body;  g  is  the  acceleration  of  gravity;  9, 
and  ^  are  the  Eulerian  Angles  for  a  ZYX  system  and  Qqj  (fioj  and  f0  are 


the  steady  flight  values, 

/ 


Defining  ,  JL  F  y  ,  and  F^y,  as  the  sum  of  the 

aerodynamic  and  thrust  forces.  Equations  3  can  be  written  as 


(9) 


£?P  = 

£Fy  -- 

i  ?L  - 

<r 


/»i-s C Ui~  (pW~~  RV)  "  P'fp' 

•bRi/ -  Pw) 

S>'nJ\  V/  ~b  PV  ~  (pU )  f 
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Each  of  the  1  'ntianeous  velocity  components  can  be  written  as  the 


sum  of  a  steady  state  value  and  a  change  in  velocity. 


U  -  Uo  + 

V  -  Vo  +  / 1/" 

W  -  V/o  tuT 

p  - 

9 r  Qotj* 

R  ~  Pofst' 

-4~-  -  ^  .^tr. 

elf 


In  order  to  linearize  the  force  and  moment  equations  it  is  assumed 

that 

1)  the  velocity  changes  are  small  and  that  the  products  and  squares 
of  the  velocity  changes  are  negligible  compared  to  the  velocity  changes 
themselves 

2)  the  disturbance  angles  are  small  and  the  sine  of  the  angles 
can  be  set  equal  to  the  angles  and  the  cosine  of  the  angles  set  equal  to 
one.  The  product  of  the  angles  is  zero. 

Equations  7  and  9  can  now  be  written  as: 


(10) 


(u) 


+  QoWo  +  Wctp t  Qoms' ~ RoVo*  Ro^- Vo/i-  "t 

Qo  )  ~ + p  (oc<u  Qp  <p*)el 

-  zw-l/lE  f  UcRo  t~  Uo si*  "t  floats  —  fo  Wo  ~  Po  -  V/of& 

OoWy  (ur^  Qo  ^  fa)  ice*'  Q0  </%)$] 

Pp  ~  srX-Mr'  -f  foVc  +  Fo/tr'  +  "  QoUo  "  i~ 

p'  Qo)  Q  ^ p-  (  Ca<u  Qo  ^y(^a)(p'^  kc<ls 

£i  ~ -tiQoRoi'Qo^'b  ~Tyy~] 
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x  M  ~  &  ]  j ')  \  0«  Ro  t*  t  ■" -L  A* J 

x  t\f  ~  y^  j  a  h  0*ci  tOo  *}  &>  y<*  f  ~  I  s/'/y'\  * 

/?  f)  £  t  /  /  ' 

The  equations  which  relate  angular  velocities  about  the  x,  y,  z  axes 
to  the  Eulerian  Angles  ©y  <P,  ^  are 

pc  q-j^e 

0  * 

Cp-  Q  Cel*-  (ft-h  y  xfc^  ft)  Co* >  & 

R  z  Qjtr&-  ft  CorZs  Q  —  ft)  , 

When  linearized  the  equations  become 


(12) 


P--  (f: 

9=  © 

r-  9 


The  steady  flight  condition  is  defined  as  all  velocity  components  equal 
to  zero  except  Up  and  W#  „  Using  this  definition  and  Equations  12,  Equa¬ 
tions  10  and  11  become 


(13) 


~  />«!-  f-  Y/op'  Qp  T Jf(c*rO  Qeft  ©J 

i  Fy  -  >«  Vat"  ■+  UoA^  ~ 

£  Fy  =  s».L  &0)O  -a  c^-doj 


(14) 

£l  = 

£  M  -  &Jyy 

In'  - 

The  aerodynamic  forces  and  moments  can  be  represented  by  a  Taylor 
series  expansion.  Because  of  the  planes  of  symmetry  and  the  steady  flight 
conditions  £.  ' <y,  ,  and  M  are  functions  of  only  q,  u,  and  w  of  the 

velocity  components  and  y  L>  and  N  are  functions  of  only  p,  v,  and  r. 


.hi  is  provt  "  i  1 «.  vce  CV  The  other  variables  affecting  the  forces 
and  i  ts  /u  ‘  the.  angle  of  thrust  deflection,  the  airframe,  flexibility^  and 
the  changing  mass  of  the  body*  Since  the  changes  are  assumed  small  the 
second  derivative  and  above  for  the  velocity  and  bending  terms  are  consider¬ 
ed  negligible  and  the  first  derivative  and  above  is  considered  negligible  for 
the  change  in  mass.  The  resulting  equations  are  of  the  form 


(15) 


where  d^  ,  ,  and  d^  are  the  first,  second,  and  third  bending  modes; 


Sr  is  the  angle  of  thrust  deflection;  and  AX  >  At  OX  ?  and  Ai-^  are 


the  changes  due  to  the  change  in  the  mass  of  the  fuels. 

The  thrust  axis  is  assumed  to  pass  through  the  center  of  gravity  and 
hence  does  not  contribute  to  the  moment  equations.  The  thrust  force  equa¬ 
tions  are 


is  the  angle  between  the  X  axis  and  thrust  line. 


T-To  -t  AT" 


and 


Therefore 


(16) 
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To  determine  t he  equations  of  steady  flight  substitute  Equations  16 
and  15  into  Equations  13  and  14t  substitute  the  values  of  the  steady  flight 
condition,  and  set  the  velocity  changes  equal  to  2ero.  The  resulting  equa¬ 
tions  are: 


(17) 


h Ta" 

s  r  O 

-~Tc  ■At+iS-  (F  Qr&'Ot  r  c> 

Lo  -O 


Mo  -  £ 

Vtf  -  (0  , 


When  stability  axes  are  used  as  the  reference  system  and  quasi-steady 
flow  conditions  are  assumed,  all  terms  containing  disappear  and  all 

aerodynamic  partial  derivatives  with  respect  to  rates  of  change  of  velocit¬ 
ies  and  with  respect  to  rate  of  change  of  angle  of  thrust  deflection  are 
eliminated.  For  the  Two  Degree  of  Freedom  Short  Period  Mode  Approximation 
MS  is  set  equal  to  zero  and  only  the  equations  F/^  and  £.M  are 
involved.  Substituting  Equations  16  and  15  into  Equations  13  and  14  and 
subtracting  Equations  17  gives  the  following  equations  of  motion  for  the 
disturbed  body: 

Jtr-  Ucfr'  t^'(^y,n^Qo)Q  -  ff^Sr  &T  +■  Fd,J/  +* 

F^l  d)  F  F  (Hz  -f  F$x  dx  4-  Fd'i  c/j  +* 

Fd'iJls  +FaxJ\i  +  FaioxAuX  +  FAlhAlH 


(IS)  #  • 

®  r  McStSt  fMd,  c£,  ft  ft  Jll  -h  Sjl  ~ft 

ft*,  ft  M/2  /i  4- Max  AT -h  MmoxALoX  -ft 

MALh  ALH 


where 


J- 


and  M 

/J 


TV7  ■j/' 
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ibc  an;  ■  t  I  1.6  d(  fined  as  the  angle  between  the  wing 

chord  line  a-.u  t  hi  r\  i<U.  '  Vt-  wind, 

yc--  ..  -  4^ 

X  ^  /\c<-  iy  -  Uv-tM'  Uo  3/tscf,  Ms-O 

(19)  Tht-ref.-  'c  lv~  Z  UoO(  • 

Equations  18v.  »  ;h  quations  12  and  19  substituted  become: 


e  =(m  US' U&  )c<  t  (aI^)  &  +(Mc {i)<$i  f  CmJ,  )<A  +(tAJj)  d* 
+  (M\M  )<^a  f  (mm)  c$3  f  (MM)  $3  -r  (m&x)AI 
+  (aj  tax  +(maih)ALH. 

The  coefficients  for  any  instant  are  constants.  Equations  20  agree  in 
form  with  the  given  equations  for  and  &  # 

The  Transfer  Functions  which  were  obtained  from  the  Aerodynamic  Equa¬ 
tions  are  sr  follows: 

1.  Aerodynamic  Condition  1; 


M-  <So 

QP  - 

*$T 

/.  ■+*.  0  6A0 

&  75- 

Qp  z 

a  •  (ksjris4'±  •  o  c>  y ) 

<£r 

i  II  y%)(r& ~t>  oAoH 9 

M-h!.9 

0P  - 

S’.oZ  (oi-.  037 H) 

<$r 

C<3-''t  Olo^J C&  — >  -t.lH  6^ 

M  -  £oo 

©p  - 

<St 

6,-,  00  ?4)C*  -/  ^i>V  &H3&) 
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j(  ’ynctiuc 


y&~c2Q 


&P 

ir 


n  2. 


SLLLL±±^M  ]C±±iQ^1  ±!&L). 

o  t>n)  (■<>  -t,  c  S’*}  t” ^  £T.  I  £> ) 


Br  -  ^17£i±±1o^£}S±^lI^A 111 

C^oU)U->  77)&+'  W) W^m>3  t#  &  A) 


X'z/o).  9 


Gp 

&T 


9  £  i  .  t>32) -t.  gjj 7  ?g,_  6, 0&2 _ 

C) )  7¥)  ('d'f',  Oft!  8.  I7y) 


zoo 


^L=  JL£2d.A±^oislL±±^3 1  zt*  a  r.i) _ 

U-,067)  Cz- /'AH)  fat  /< 3 (*,+■,  a7,  7  ) 


3.  Aerodynamic  Cord  it  ion  3; 


/:  30  JSf  =  /6.q2Tc *  + 

f,  oft  tp  fc  ttf) 


j£  &R  -  1 9.  /  (<ch- «  V-,  a  <$  7  1 1 z^L  QJlL 4 ).  _ _ 

^  e^T  0*i  ~ ,  o<A 3^)  ~i  7 7^J V-,  8 >  I ®3 Bi HS"3  HI- 


/  -  /c/,  9  -  siL  zLzlzLZ.  >£kji&dl  2££.  -  92JL _ 

&r  (<!y-,c  I )(<$,  -.  &83j  Ca  i.  7  0  87-^8".  !Zj(&-b  Ul<63) 
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c2oo 


&£  z  VU+- 6 &)(<*  1 1  il  1  A r.  73,3) 
^T  ^  ',to?)(*-/.MK-*+/>  ■>)(*+'  m  ty  73, 1 ) 


Aerodynamic  Condition  4; 


ri/^  &P  -  to.  78  Cat  •0^)(^i'.0lS'tXl,  H%)Uf.olc/'i^-/' 69)C4-+<o£&tiy  J>63) 


( d±A  %  -.?  !L  7  3)  (o  z  lh  sll  t^h  Ull&d.  £h  &J. 

(<0-+ ,018'i‘J  /, 9£>)  i^f<OlS‘ t^.' /  Si)  (a. ,  IGH) 


-Zz  IS  ®p  -  / ?. I  U+. c?)(^t,oo7 -*V /  yQQp r> QM  tf  J.  ,2 ?)Lt, QZ 7 ?? X, C.Q 

~%r  (^tAsJ  U-.ol)U- ,  ta)  U^.  ysi £  tjlhir) 

(  ,j /•■  c  9H  -  X-9. 37^( ~ i6<  A  —  X  rOAi  ~?S^) C&s+lGi  C?  t  *y  A2  < 

(Zt\  06  ± jf/o)(^+.  }o3  tJ 


A-z/qI.J  £Lp  z  & 1  rttH)U,  ~%3 jj^MxJ£s2^Lt C/-  —  tt^2 - 

St*  -idOCd-  ~ '  ^3  » H  £\%5£j 

f^r .  ±>  H.Al  )L  +.CG  r ft,  >  )%)(&+.  032. 

(^t,03C  Th  A  H)  uT°*Tfy  8.  Is)  U  +  ,03H±  p  3 ,  HZ) 


y~  -  ,  (~9p  -  sLl-SjSI L&3Li  'l *A7 ~jy A ■  3a 

*5r  (*  - , ca  v ;  ^  - .  /r.ri^  /.  7^)^+.  7a¥±^  Try?/ 

X^Ju  HSS  Zj,  V  (>.  7^)U  ti  *£&  73,  l) _ _ 

C *+,A97  t/c2Z  %)C++^i  r  f^-  HLl)(or,7HSt^3,S) 
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■ 


SECOND  BENDING  MODE 


Computer  Simulation  Diagram. 
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Figure  B.l.  Continued. 


>C— o'  -6^-0 
« 


98 


Figure  B.l.  Continued. 


Figure  B.l.  Continued. 
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Figure  B.2(a).  Frequency  Response  Amplitude  vs  for  Dither  Frequency 

of  50  ana  Aerodynamic  Condition  1.  FHt 
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ANAI  YXI^AL  VZHrSTrAITOK  OF  xAL  EiVIl-GCAi  RECTI’ 

The  results  of  the  variation  of  damping  ratios  with  dither  frequency 
obtained  by  empir:  ?A  methods  in  Chapter  4  will  be  analyzed  by  analytical 
methods  in  this  Appendix* 

For  the  analysis  the  root  loci  for  the  flight  cases  of  a  particular 
aerodynamic  cond  tion  will  be  assumed  to  have  the  same  shape  and  system 
gain  values.  The  transfer  functions  of  the  vehicle  dynamics  for  a  parti¬ 
cular  aerodynamic  condition  do  not  vary  to  the  extent  that  the  use  of  the 
assumption  to  determine  general  trends  would  be  prevented*  A  particular 
system  gain*  K,  on  the  root  loci  for  the  three  flight  cases  of  a  condition 
xvould  determine  the  same  damping  ratio^  The  system  gain  is  the  produc 
of  the  dynamic  gain,  M £  ,  and  the  adaptive  gain,,  K£  0  The  dynamic  gain 

is  the  gain  associated  with  the  vehicle  transfer  function*  Using  condition 
3  as  an  example ,  the  dynamic  gains  for  the  flight  cases  are  approximately 
equal  to  the  following  values: 


Case  t  20 

vs  ~- 

16.05 

Case  t  7 5 

19. 1 

Case  t  101.9 

21.6 

Since  a  particular  damping  ratio  has  the  same 

va  1  ue 

of  R  for  the  three 

cases,  the  values  of  for  the  cases 

are : 

K 

K  &■  Z  -  2  0 

16.05 

k 

H/r=  7iT 

19.1 

K 

r  Jol.9 

21.6 

and  Kst  20  >  K  £  t  -  75  >  K£ 

t  - 

101.9 

for  a  particular  £ 

A  plot  of  ^  versus  £  is  of  the  following  form„ 
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ks 


The  trends  shown  in  the  <£  versus  K,£  plot  are  valid  for  any  of  the  four 
aerodynamic  conditions. 

The  closed  loop  transfer  function  for  the  system  is 

e  _  k^A  _ _ 

ceMM/j/ZD  Si&rfAL  ~  li-Kz(^  +  3.6«2,+  3.6>) 
when  *-s  the  vehicle  dynamic  transfer  function.  When  the  dither  freq- 

uency  is  large  compared  to  the  poles  and  zeros  of  the  vehicle  dynamic  trans¬ 


fer  function,  the  transfer  function  is  approximated  by 


&£ 


St 

This  approximation  is  justified  for  the  dither  range*  10  to  40  radians  per 
second,  for  conditions  1  and  2.  The  closed  loop  transfer  function  becomes 


© 


Cot* 1MAJ/D  SIGNAL 


}+  k<?  M± 


Since  the  dither  system  maintains  the  amplitude  of  the  closed  loop  response 
constant „ 


K&,  Mifi 

-  kfiM&t 

-  KsjlMsi 

1  i  k±,  Ms, 

jt-ao  1-HfoMw 

£.7S  H-K^Mss 

and  for  Aerodynamic  cuuuition  1 
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in 


kti  -  KxA  J.  Ujj Z  - 

)  -h  Ks,  i,  SCS  1-hK^J.  US'  !  +  Ksi  S.oQ> 


and  for  Aerodynamic  Condition  2 

to,  £_  IS  -  J (IulAJS:  ^  kss  9.  A 

It-kjiSjs  H-Ksj6'7S'  H-K*i?,A 

The  relationships  are  independent  of  the  dither  frequency*  and  constant 
damping  ratio  variation  would  be  expected  in  the  dither  range.  This  agrees 
with  the  empirical  results  for  conditions  1  and  2  as  are  shown  in  Figure 
4.4.  If  is  assigned  the  value  1  for  both  condition  1  and  2^  then  for 

condition  1 


K£  {  =  1,67  K£3  =  0,516 

and  for  condition  2 

!-3  KS  3  "  °'725 

For  both  conditions,  ^  ;  )  K  ^  which  is  the  trend  necessary 

for  small  damping  ratiw  variation  as  determined  from  the  cf  versus  plot. 

For  condition  3  the  dither  frequency  can  be  close  to  poles  and  zeros 
of  the  vehicle  dynamic  transfer  function.  If  the  dither  frequency  is  close 
to  a  pole  of  the  transfer  function,  then  is  very  large  and 

k&C £+3-6*+  3>  0f£>?l 


G 


and 

A 


coMwAr/o 

The  amplitude  of  the  closed  loop  response  is  independent  of  the  adaptive 
gain;  and  therefore^  I<£  could  assume  a  infinite  number  of  values  with 
an  infinite  number  of  damping  ratio  variations. 

If  the  dither  frequency  Ls  close  to  a  zero  of  Case  t  =  20,  then 


e, 


< 


< 


Bs. 


Sr^yjr  &Tj£-  I  oh  f 
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}  Ks  «■  -  }  K* 

variation.  If  the 


s'  anplifuut.  of  T  he  closed  loop  re  sp^nse*  t  =  20 

101.9  which  Is  the  trend  for  a  small  damping  ratio 
dither  frequency  is  close  to  a  zero  of  Case  t  ~  101,9, 


then 

Op  v  Bf  \  &£  , 

'  S>t£z7£‘  '  Sr^ilal,1) 

In  order  to  have  the  same  amplitude  of  the  closed  loop  response  Ks  t  =  20< 
Ks  r  =  75  ~  101.9  which  is  the  trend  which  would  result  in  a 

large  damping  ratio  variation.  If  the  dither  frequency  is  close  to  a  zero 
of  Case  t  =  75,  the  magnitude  of  the  damping  ratio  variation  would  vary 
between  the  magnitudes  when  the  frequency  is  close  to  a  zero  of  Case  t  ~  20 
and  when  the  frequency  is  close  to  a  zero  of  Case  t  =  101,9* 

The  damping  ratio  variation  curve  for  condition  3  should,  therefore, 
have  an  undetermined  damping  ratio  variation  at  the  pole  at  approximately 
40  radians,  have  a  small  value  for  the  variation  at  the  zero  of  case  t  -  20 
at  approximately  24  radians;  and  have  a  large  value  for  the  variation  at 
the  zero  of  case  t  =  101.9  at  approximately  27  radians.  This  agrees  with 
the  condition  3  curve  shown  in  Figure  4*4, 

In  order  to  determine  the  dither  frequency  that  the  minimum  variation 
in  damping  ratio  occurs,  the  ideal  expression  for  no  variation  which  is 
of  the  form 


Jr  3r 
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sr 
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wbe  n_ 


n  ~  ^  s  1 

3  =  KJ  t 
C  =  KJ  t 


75 


J  f  30 

go 

sJrfJO 


101.9  ^ 

•d^f  Jo 


and  for  a  damping  ratio  of  0,7 
K  £  r  -  20  -  0.86 
Kf  t  =  75  =  0,26 

<y 

t  -  101*9  =  0 *  175 

is  evaluated  for  various  dither  frequencies  close  to  the  frequency  of 
the  zero  for  case  t  =  20,  and  the  smallest  amplitude  variation  between  the 
three  cases  should  indicate  the  optimum  frequency*  A  dither  frequency  of 
23  radians  per  second  resulted  in  the  smallest  variation  in  the  amplitude. 

For  aerodynamic  condition  4  there  are  no  zeros  or  poles  of  the  vehicle 
dynamic  transfer  function  which  are  close  to  the  imaginary  axis  in  the 
dither  range.  Varying  the  dither  frequency  does  not  result  in  wide  varia¬ 
tions  in  the  values  of  the  vehicle  transfer  functions*  and  the  closed  loop 
transfer  function  can  be  approximated  by 


e 

ComMA/ZO  signal 


M± 

)+  k 


A  constant  damping  ratio  variation  would  be  expected  and  agrees  with  the 
curve  for  condition  4  in  Figure  4.4.  Since  the  values  of  the  vehicle 
transfer  functions  can  not  be  varied  appreciably,  the  variation  in  the 
adaptive  gain  required  for  a  constant  damping  ratio  can  not  be  realized 
and  a  large  variation  in  the  damping  ratio  results* 
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